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ABSTRACT 

We present an analysis of the maximum light, near ultraviolet (NUV; 2900 < A < 5500 A) 
spectra of 32 low redshift (0.001 < z < 0.08) Type la supernovae (SNe la), obtained with 
the Hubble Space Telescope (HST) using the Space Telescope Imaging Spectrograph. We 
combine this spectroscopic sample with high-quality gri light curves obtained with robotic 
telescopes to measure SN la photometric parameters, such as stretch (light curve width), op- 
tical colour, and brightness (Hubble residual). By comparing our new data to a comparable 
sample of SNe la at intermediate redshift (0.4 < z < 0.9), we detect modest spectral evolution 
(3-cr), in the sense that our mean low redshift NUV spectrum has a depressed flux compared 
to its intermediate redshift counterpart. We also see a strongly increased dispersion about the 
mean with decreasing wavelength, confirming the results of earlier surveys. We show these 
trends are consistent with changes in metallicity as predicted by contemporary SN la spectral 
models. We also examine the properties of various NUV spectral diagnostics in the individual 
SN spectra. We find a general correlation between SN stretch and the velocity (or position) 
of many NUV spectral features. In particular, we observe that higher stretch SNe have larger 
Can H&K velocities, which also correlate with host galaxy stellar mass. This latter trend is 
probably driven by the well-established correlation between stretch and host galaxy stellar 
mass. We find no significant trends between UV spectral features and optical colour. Mean 
spectra constructed according to whether the SN has a positive or negative Hubble residual 
show very little difference at NUV wavelengths, indicating that the NUV evolution and varia- 
tion we identify does not directly correlate with Hubble diagram residuals. Our work confirms 
and strengthens earlier conclusions regarding the complex behaviour of SNe la in the NUV 
spectral region, but suggests the correlations we find are more useful in constraining progeni- 
tor models rather than improving the use of SNe la as cosmological probes. 

Key words: distance scale - supernovae: general - galaxies: general 
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1 INTRODUCTION 

The use of Type la supemovae (SNe la) as cosmological probes is 
now well established; significant advancements i n constraining the 
cosmological parameters using distant SNe la (Rie ss et al, I l2007l : 
Kessler et al. 20091; ISullivan et alj|20l l|) have been made since the 
initial di scovery of the accelerating Universe at th e end of the last 
century jRiess et alj[l998l ; IPerlmutter etaill 19991) . SNe la are be- 
lieved to result from the thermonuclear explosion of an accreting 
carbon-oxygen (CO) white dwarf in a binary system. Although the 
nature of the p rimary exploding star has recently been confirmed as 
a white dwarf (Nugen Fet al.l201 lllBloom et al.l2012T) . the nature of 
the secondary, mass-donating companion star remains a mystery. 
Generally, two possibilities are considered - double-degenerate 
systems (with two white dwarfs) and single degenerate scenarios (a 
white dwarf plus a giant, sub-giant or main sequence star). Recent 
results have indicated a diversity in the form of this companion star 
from SN event to SN e vent, with evidence that some SNe la have 
red giant companions jPildav et alj|2012 | ). evidence tha t directly 
exclu d es red-giants in other cases ([Nugent et all l201lt lli et alj 
201 ll ; IChomiuk et alj I2OI2I; |Horesh et al] |2012| ; iMargutti et"aLl 
iRuss 



2012 ; [Russell &Immlej 120121)? eviden ce that favours single de- 



generate systems (Sternberg et al. 201 1), and evidence that favours 
double-degenerates dSchaefer & Pagnottall2012l) . If SNe la do in- 
deed result from two (or more) progenitor channels, then under- 
standing their differences and their dependence on host galaxy 
properties (or stellar populations) is of utmost importance for their 
use as distance indicators. 

The metallicity or composition of the exploding white dwarf 
is also likely to affect the photometric properti es of SNe la 
dTimmes. Brown & Truranll2003t iKasen et alj|2009l) . an effect that 
may have been detected observationally: Hubble residuals correlate 
with the stellar mass of the SN la host galaxies jKellv et alj|2009l ; 
ISullivan et alj20ld ; lLampeitl et alj201ol) . which can be interpreted 
as a crude proxy for metallicity. Gas-phase metallicity measure- 
ments have also been found to correlate with Hubble residuals 
dD' Andrea et alj|201 lb . SN la spectroscopy is a potential probe of 
this astrophysics, and the rest-frame ultraviolet (UV) in particular 
is expected t o show strong signatures of metallicity/compositional 
effects (e.g., Hoflich, Wheeler and Thielemann 1998; iLentz et al] 
l200d;ISaueVetai]|2008l ; IWalker et al.ll2012l) . The NUV spectral re- 
gion is dominated by numerous overlapping lines, that make identi- 
fying individual features difficult and results in heavy line blanket- 
ing. The UV flux of SNe la at early times is formed mainly through 
the process of 'reverse-fluorescence' , the scattering of photon s 
from longer to shorter wavelengths dLucvl [T999t lMazzalil l2000). 
The UV photons emitted through this process are predominantly 
found to ori ginate from a thin layer of material at high-velocity 
dMazzalil2000l) . Therefore, the NUV spectral region is a good probe 
of the outer layers of the SN ejecta. 

Generally, spectral studies have taken two different ap- 
proaches: either a comparison of mean spectra, testing for evolu- 
tionary effects in composition or velocity profile with redshift, or 
measurements of individual spectral features, testing for relation- 
ships between these spectral features and SN photometric proper- 
ties. Some spectral comparisons between low- and high-redshift (z) 
samples have suggested that there may be some evolution with red- 
shift i n the mean NUV spec t rum (Foley et a l. 2008a; B alland et al.1 
120091 ; ISullivan et all 120091 ; iFolev et al] l2012j ; lEllis et al] |2008|. 
hereafter, E08). However, due to the difficulties in obtaining UV 
spectra of low-z SNe la, studies of the metallicity-influenced UV 
region have been limited by a small number of spectra of just a 



handful of SNe la at low redshift, and have not been able to probe 
either the full range of "normal" SN la properties or host galaxy 
properties. Observations at higher redshift - where the rest-frame 
near-UV (NUV) is redshifted into the optical and observations are 
dominated by UV flux - have also shown an increased scatter 
at NUV wavelengths compared to the optical (E08). The limited 
Hubble Space Telescope (HST) and SWIFT spectrosc opy of low-z 
events that is available supports the higher - z studies dBufano et al] 
l2009l ; ICooke et al]|201ll : IWang et all2012l) . as does Swift UV pho- 
tometry which illustrates an incre ased dispersion in UV colours 
imply ing larger spectral variations dMilne et al.boiot iBrown et al] 
2010). One interpretation is that the larger variations at these wave- 
lengths are due to differing progenitor metallicity and composition; 
an increased dispersion is in broad agreement w ith the results of 
SN la UV spectral modelling dWalker et al]|2012l) . 

Studies of particular SN la spectral features have usually con- 
centrated on improving their use as distance indicators, either by 
measuring spectral feature equivalent widths, velocities, or spec- 
tral flux ratios. Various empirical trends of differing statistical sig- 
nificance have been claimed be tween SN light-curve width and 
spect r al feature equivalent width dBronder et al.l20 08; Walk er"et al] 
1201 1; Chotard et al. 2011) , between SN luminosity and spectral 
flux ratios ([Foley et all l2008bl ; iBailev et al] 120091; iBlondin et all 
l201ll ; Isilverman et al]|2012al), and between SN colour and spec- 
tral velocities |Wang et al]|2009l : [Foley. Sanders & Kirshnerll201ll ; 
Blond in et al 1 l20 1 2l ; lFolevll20 1 2h . As of yet, none of these spectral 
trends have significantly out-performed the standard photometric 
indices used in constructing Hubble diagrams. Trends between fea- 
ture velocity and host galaxy stellar mass have also been identified 
dFoley|20"l2l) that may hint at a connection b etween spectral proper- 
ties and the underlying progenitor systems. Sternberg et al] J20T1I) 
recently provided evidence of systematically blue-shifted absorp- 
tion components in SNe la, likely related to pre-explosion outflow, 
favouring the single-degenerate proge nitor scenario in s piral galax- 
ies. This work has been extended by IFolev et al] d2012bl) to show 
that the SNe with blue-shifted circumstellar absorption features, on 
average, display higher ejecta velocities and redder maximum light 
colours than normal SNe la. 

The above results and correlations are very important, both 
in our understanding of the physical diversity of the SN la phe- 
nomenon as well as the potential for biases in the use of SNe la as 
cosmological probes. However, the local UV samples have been, to 
date, modest in size. In this paper, we rectify this shortcoming. 

Our goals are several. Using a new, unbiased, low-z sample 
we aim to verify whether there is evolution with redshift in the 
NUV spectra, as well as investigate the claimed correlations be- 
tween NUV spectral, light curve and host galaxy properties, and 
their effect on the Hubble residuals. To do this, we use the maxi- 
mum light spectra of 32 low redshift SNe la (0.001 < z < 0.08) 
at wavelengths down to ==2900A (NUV) obtained using the HST. 
Our new sample e xpands the initial study of 12 events presented in 
ICookeetaljfcOllI) , and now also includes additional multi-colour 
light curve data, which are used to colour-correct the spectra. All of 
our NUV spectra were obtained with the Space Telescope Imaging 
Spectrograph (STIS) on HST; 28 were obtained during HST Cycle 
17 (GO 11721, PI: Ellis), and 4 during Cycle 18 (GO 12298, PI: 
Ellis). Although data down to ~3200 A can be obtained from the 
ground, these data are notoriously difficult to calibrate near the at- 
mospheric cut-off. The two key advantages of HST are an accurate 
relative flux calibration over the entire wavelength range, and min- 
imal host galaxy contamination of the spectra due to the narrow slit 
employed. 
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We primari ly select our SNe la from the P alomar Transient 
Factory (PTF. iRau etafl 120091 ; IL"aw et alj 12009). a low-z rolling 
search that selects SNe in a similar manner as those found in high-z 
surveys, without any bias towards a particular host type. We use 
these data to compare to the NUV spectra of 36 intermediate- 
redshift (0.4 < z < 0.9) SNe la from the Supernova Legacy Survey 
(E08) to look for evolution in the SN la properties with redshift and 
compare the dispersion of their spectra, for the first time using an 
unbiased low-z sample. This homogenous low-z sample can also be 
split into mean spectra based on their light curve and host galaxy 
properties. The positions (velocities) of the NUV spectral features 
identified in the individual spectra can also be measured qualita- 
tively and compared with their stretches, colours, Hubble residuals 
and host galaxy properties. 

The plan of the paper follows. § [2] introduces the data used 
in this paper, including the HST spectra, the ground-based optical 
light curves, and the SN la light curve fitting technique. This sec- 
tion also gives information on the host galaxies and the method 
for the calculation of Hubble residuals. § [3] details the results of 
the mean spectrum comparisons, including those split by redshift, 
as well as SN and host galaxy properties. §|4]describes the results 
of wavelength and velocity measurements of spectral features and 
how they relate to observed SN quantities and Hubble residuals. 
The results of the previous two sections will be discussed in § [5] 
including the implications for the use of NUV spectral features as 
further calibrators of the Hubble diagram, how these correlations 
relate to possible progenitor scenarios, and the cause of the ob- 
served evolution with redshift and dispersion at NUV wavelengths. 
Throughout this paper, where necessary we assume a Hubble con- 
stant, H = 70kms~' Mpc -1 . 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sample selection 

The bulk of our NUV spectra of SNe la come from a 'non- 
disruptive' Target of Opportunity (ToO) programme (GO 11721, 
PI: Ellis). The scheduling constraints of HST meant that the SN 
la targets had to be submitted for scheduling ~7— 16 d before the 
required m aximum-light spe ctrum (for further details of our tech- 
niques, see ICookeetal.|[201lh . This requires very early detection 
and spectroscopic classification of potential candiates, which were 
discovered by the PTF, and classified by a mixture of scheduled and 
ToO spectroscopic time on a large number of optical telescopes, de- 
tailed below. 

PTF is a wide area, rolling search optical survey using the 
Samuel O schin 48-in Tele s cope (P48) locat ed at the Palomar Ob- 
servatory teauetal .1120091 11" aw et al. 2009). Variable cadences in 
the range of minutes to ~5 d are routin ely used to identify S N can- 
didates very soon after explosion (e.g. iGal-Yam et al 1 l201ll) . 27 of 
the 32 low-z SNe la studied here were discovered by PTF, with 8 of 
the PTF SNe being ide ntified using the " Galaxy Zoo: Supernova" 
citizen science project dSmithetal.ll201lh . Five other SNe la were 
discovered by ot her searches: SN 2 0091e (Chilean Automatic Su- 
pernova Search: IPignata et al.l2009h. SN 2010 iu (Lick Observatory 
Supernova Search (LOSS); ICenko et al.l2010l) , SN 2010kg (LOSS; 
iNavak et al.ll2010h. SN 201 lby (Xin gming Observation Sky Sur- 



Jin et alj|20lil) and SN 201 lek jNakano et al]|201 ll) . Tabled 



vey; 

details the optical classification spectra of the sample, the redshifts 
of the SNe and their host galaxy stellar masses (M stc i lar , which will 
be discussed further in SectionlZTot, 



30 NUV spectra of 30 different SNe were obtained during 
the main cycle 17 programme, 28 of which were "normal" SNe 
la, one of which w as the unusual SN la PTFlOops (detailed in 
Magu ire~et al.ll20lTI) . and one of which was ultimately classified 
as a SN Ic (PTFlOosn). These latter two events are not studied in 
this paper. During HST Cycle 18, a follow-on HST programme (GO 
12298, P.I.: Ellis) studied the time evolution in the NUV and far- 
UV spectra of a smaller sample of SNe la from soon after explosion 
to after maximum light. This programme is complementary to the 
Cycle 17 programme and will be presented in detail in further pub- 
lications. However, we have supplemented our sample of the 28 
"normal" cycle 17 SNe la maximum light spectra with four maxi- 
mum light spectra from this progr amme: PTFlOygu ( Hachinger et 
al. in prep.), PTF1 lkly/SN 201 lfe dNugent et alj201lb . SN 201 lek 
and SN 201 lby. This brings our sample to a total of 32 maximum- 
light NUV spectra, by far the largest NUV spectral sample of low-z 
SNe la to date. 

The light curve derived properties (date of maximum, B-band 
magnitude at maximum, stretch and B - V colour at maximum) of 
the sample are described in Table [2] We will compare our sample 
to that of E08, which contains a large sample of intermediate-z SN 
la UV spectra that were obtained using Keck+LRIS. Their sample 
of SNe la was d iscovered by the Supernova Legacy Survey (SNLS; 
lGuvetal.| [2010) and are in the redshift range of ~0.4-0.9, with a 
mean redshift of 0.6. We extend the intermediate-z sample of E08 
to include nine further SNLS SNe la spectra obtained with the same 
Keck setup (although not part of the E08 sample) and reduced using 
the same method as described in E08. Throughout this paper, refer- 
ences to the intermediate-z E08 sample will include these nine ad- 
ditional SN spectra. A comparison of the effective phase (the phase 
of the NUV spectrum divided by the SN stretch), stretch, B - V 
colour at maximum, and M ste uar distributions of the low-z sample 
to that of the intermediate sample of E08 is shown in Fig. [T] along 
with the full three-year SNLS SN la sample (SNLS3: Tguv et all 
l20ld: IConlevetalj201ll: ISullivanetal]201ll) . 

Using a Kolmogorov-Smirnov (K-S) test, we find that there 
is a very low probability of our low-z HST sample and the 
intermediate-z sample of E08 having their effective phase, stretch, 
B - V colour at maximum, and M stc n ar distributions drawn from 
different parent populations (using a phase range of -1.0-4.5 d, the 
choice of which will be discussed in Section [3~0 . This suggests that 
the chosen low- and intermediate-z samples are well-matched and 
any evolution in their spectral properties should not be due to dif- 
ferences in the phase and light curve properties of the samples be- 
ing studied. Luminosity biases that could occur because of brighter 
SNe la being discovered preferentially at intermediate-z would be 
manifested as a stretch or colour bias, which is not seen. 



2.2 UV spectroscopy 

The HST spectra were obtained using STIS and the G430L grism, 
giving a wavelength coverage of 2900-5700A with a dispersion of 
2.73 A per pixel and a plate scale of 0.051" per pixel. Fig.[2]shows 
the NUV spectra which are discussed in this pap er. All the spectra 
are pu blicly released via the WISeREP portajj] dYaron & Gal- Yam] 
l2012h . 

The spectr a were reduced in a manner broadly following that 
of ICooke et all feoill) . The spectra were downloaded from the 
HST archive using the on-the-fly reprocessing (OTFR) pipeline. 



http://www.weizmann.ac.il/astrophysics/wiserep/; 
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Figure 1. Comparisons of the SN la properties for the different samples used in t his paper. In the t op left panel, the effective phase (the SN phase divided by 
stretch) distribution of the sample in this paper (gray), along with the sample of lEllis et al] 1200 8l) (E08, red) is shown. The vertical, dashed lines mark the 
positions of the phase cuts at -1.0 and +4.5 d that are used to choose the SNe la that enter the mean spectra. The top right and bottom left panels show the 
stretch and B ~ V colour distributions, respectively, for the phase-selected low-z sample (gray), the phase-selected sample from E08 (red) and that of the full 
SNLS3 SN la sample from lGuv et al . (2010) (G10, blue). These two plots have had their distributions normalised to one for a clearer comparison. The lower 
right panel shows the host galaxy stellar mass distributions, (M stc ih u ■), respectively for the phase-selected low-z sample (gray) and the sample of E08 (red), as 
calculated from fits to the host galaxy magnitudes. If the host galaxy was not detected (PTFlOufj), a value of M stc u al . of 10 6 M G was assigned. 



This gives fully calibrated ID spectra ('sxl' files), where the re- 
duction and extraction is optimised for point sources. The OTFR 
pipeline uses the latest calibration files and data parameters to per- 
form initial 2D image reduction such as image trimming, bias and 
dark current subtraction, cosmic-ray rejection via CRSPLIT and 
fiat-fielding. It then performs wavelength and flux calibrations, in- 
cluding correcting the flux for imperfect charge transfer efficien- 
cies along the c hip. We applied fur ther cosmic ray removal us- 
ing LACOSMIC (Van Dok kumll200ll) and corrected the spectra for 
Milky Way ex t inctio n using an R v value of 3.1, the dust maps of 
ISchleg el et al. ( 1998), and the Milky Way dust extinction relation 
of lCardelli. Clayton & Mathisl dl98St CCM). Most of the SNe la in 
the sample had low Milky Way extinction values of E(B-V) < 0.2, 
although two, SN 201 lek and SN 2010ju had values of 0.35 and 
0.42 respectively. SN 201 lek had an unusually red UV minus opti- 
cal colour and the HST spectrum of SN 2010ju is outside the phase 
range discussed in the rest of the paper. Both have been excluded 
from further analysis. 

We assess the possibility of host galaxy contamination in our 
HST/STIS spectra using the acquisition images taken with each 



spectrum. These use the STIS/CCD long-pass filter and have ex- 
posure times ranging from 5-35 s, depending on the apparent mag- 
nitude of the SN. In all cases, the SN is by far the brightest object 
in the image, and only 7 of the images show any visual trace of 
the host galaxy. For each SN, we estimate the possible contamina- 
tion by placing circular apertures of 0.2" diameter (the slit width, 
approximately 4 STIS/CCD pixels) at the position of the SN, and 
also along the circumference of a circle of radius 0.6" centered on 
the SN position. These latter apertures measure the host galaxy flux 
near to the SN position but are outside the extent of the STIS/CCD 
point spread function (PSF). We take the brightest of these host 
galaxy apertures as conservatively representing the flux attributable 
to the host galaxy. We then calculate the percentage excess in the 
brightness of the SN over its host. For 30 SNe (94 per cent of the 
sample), the SN light through the slit is more than 100 times the 
estimated brightness of the host through the slit. Even for the worst 
case SN (PTFlOxyt), the SN flux through the slit is -30 times that 
of the host galaxy flux. Therefore, we consider the host galaxy con- 
tamination to be negligible in our analysis. 

Spectral line identifications were made for all the HST spec- 
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Table 1. Log of optical classification and HST spectra, along with the SN redshifts and host galaxy log(M stc u al -) values. 



SN name 


R.A. 


Decl. 


Class. Telescope/ 


MJD 


MJD HST 


CMB redshift (z) c 


Log(M sKllar ) 




(J2000) 


(J2000) 


CBET 


Classification 


spectrum 




(M G ) 


PTF09dlc 


21:46:30.10 


+06:25:09.2 


P200+DBSP 


55063.5 


55076.7 


0.0666±O.00O5 


9.01 


PTF09dnl 


17:23:41.79 


+30:29:49.6 


P200+DBSP 


55063.5 


55076.1 


0.019±0.001 


7.84 


PTF09dnp 


15:19:24.37 


+49:29:56.3 


P200+DBSP 


55063.2 


55076.6 


0.037559±0.0000l' 


11.10 


PTF09fox 


23:20:47.98 


+32:30:08.2 


Keck+DEIMOS 


55123 


55134.1 


0.0707±0.0005 


10.39 


PTF09foz 


00:42:11.73 


-09:52:52.7 


P200+DBSP 


55126 


55134.2 


0.05331±0.00002 


10.94 


PTFlObjs 


13:01:11.23 


+53:48:57.3 


WHT+ISIS 


55250.2 


55263.7 


0.0303±0.0003 


6.85 


PTFlOfps* 


13:29:25.06 


+ 11:47:46.5 


P200+DBSP 


55303 


55319.5 


0.0225 13±0.000009' 


10.90 


PTFlOhdv 


12:07:45.37 


+41:29:27.9 


KECK+LRIS 


55331 


55347.6 


0.0542±0.0003 


7.50 


PTFlOhmv* 


12:11:32.99 


+47:16:29.8 


P200+DBSP 


55337 


55354.3 


0.033±0.001 


8.34 


PTFlOicb 


12:54:49.21 


+58:52:54.7 


Gemini N+GMOS 


55349 


55361.4 


0.0088±0.0003 


9.66 


PTFlOmwb* 


17:17:49.97 


+40:52:52.1 


Gemini N+GMOS 


55367 


55390.4 


0.0312±0.0005 


9.36 


PTFlOndc" 


17:19:50.18 


+28:41:57.5 


P200+DBSP 


55384 


55396.5 


0.0817±0.0005 


9.05 


PTFlOnlg 


16:50:34.48 


+60:16:35.0 


Keck+LRIS 


55384.4 


55397.3 


0.056±0.005 


10.11 


PTFlOpdf 


13:24:30.98 


+42:05:20.3 


P200+DBSP 


55396 


55410.1 


0.0764±0.0005 


11.08 


PTFlOqjl 


16:39:59.33 


+12:06:25.5 


WHT+ISIS 


55411.0 


55424.4 


0.0579±0.0005 


9.10 


PTFlOqjq 


17:07:12.39 


+35:30:35.3 


WHT+ISIS 


55411.0 


55424.4 


0.0288±0.0005 


10.14 


PTFlOqyx* 


02:27:12.06 


-04:32:04.8 


Lick+KAST 


55419.5 


55431.8 


0.065±0.005 2 


6.00 


PTFlOtce"- 6 


23:19:10.36 


+09:11:54.2 


Gemini N+GMOS 


55433 


55445.7 


0.0397 16±0.000017 3 


10.60 


PTFlOuff-' 


02:25:39.13 


+24:45:53.2 


P200+DBSP 


55444 


55459.1 


O.076±0.005 2 


6.00 


PTFlOwnm 


00:22:03.61 


+27:02:26.2 


VLT+XSHOOTER 


55467 


55480.6 


0.0645±0.0001 


10.52 


PTFlOwoF 


23:32:41.84 


+ 15:21:31.7 


Gemini N+GMOS 


55468 


55480.1 


0.0514±0.0005 


10.18 


PTFlOxyt" 


23:19:02.42 


+ 13:47:26.8 


KPNO+RC Spec 


55480.5 


55494.3 


0.0484±0.0003 


9.75 


PTFlOygu"* 


09:37:30.29 


+23:09:33.0 


Gemini N+GMOS 


55482 


55495.4 


0.0260±0.0005 


11.20 


PTFlOyux" 


23:24:13.39 


+07:13:42.7 


Gemini N+GMOS 


55487 


55502.3 


0.0559±0.0005 


11.22 


PTFlOzdk" 


02:14:07.33 


+23:37:50.2 


WHT+ACAM 


55502.0 


55515.1 


0.032±0.001 


9.70 


PTFlOacdh 


09:43:07.58 


+09:39:31.3 


P200+DBSP 


55543 


55558.1 


0.0606±0.0002 1 


10.82 


PTF1 lkly 


14:03:05.83 


+54:16:25.2 


LT+FRODOSPEC 


55797.9 


55814.4 


0.001208±0.000005 




SN20091e 


02:09:17.14 


-23:24:44.8 


CBET 2022/2025 


55151 


55165.2 


0.01703±0.000030 3 


11.04 


SN2010ju 


05:41:55.99 


+18:29:51.0 


CBET 2549/2550 


55516.3 


55529.6 


0.015348±0.000013 3 


10.93 


SN2010kg 


04:40:08.40 


+07:21:00.0 


CBET 2561 


55531 


55543.3 


0.013968±0.000013 3 




SN201 lby 


11:55:45.56 


+55:19:33.8 


CBET 2708 


55678.5 


55690.3 


0.003402+0.000005 3 


9.27 


SN2011ek 


02:25:48.89 


+18:32:00.0 


CBET 2783 


55780.2 


55792.4 


0.004206+0.000010 3 





"'Galaxy Zoo Supernovae project' discovered SNe. 

*IAU designated name, PTF10fps=SN 2010cr, PTF10hmv=SN 2010dm, PTF10mwb=SN 2()10gn, PTF10qyx=SN 2010gy=PSl-100385, PTF10tce=SN 
2010ho, PTF10ufj=SN 2010hs, PTF10ygu=SN 2010jn 

c lf not otherwise marked, z values are measured from host galaxy features . 'SPSS, template fit, 3 NED 

P200+DBSP: Palomar 200-inch telescope with the Double Spectrograph jpke & Gunnll982h. 

Keck+DEIMOS: Keck 10-m with the Deep Imaging Multi-Object Spectrograph jFaber et~ 20031) . 

WHT+ISIS : 4.2-m William Herschel Telescope with the Intermediate dis persion Spectrog raph and Imaging System. 

Keck+LRIS: Keck 10-m with the Low-Resolution Imaging Spectrometer (Oke et al. 1995). 

Gemini N+GMOS: Gemini North with the Gemini Multi-Object Spectrograph (Ho ok et alj2004l) . 

Lick+KAST: Lick Observatory Shane 3-m telescope with the KAST double sp ectrograph. 

VLT+XSHOOTER: Very Large Telescope with the XSHOOTER spectrograph iD'Odorico et al.l2006tlvernet et all20TTh . 
KPNO+RC Spec: KPNO 4-m Mayall telescope with the RC spectrograph. 
WHT+ACAM: 4.2-m William Herschel Telescope with the Auxiliary-Port Camera. 
LT+FRODOSPEC: Liverpool Telescope with the Fibre-fed Robotic Dual-beam Optical Spectrograph. 



tra in the sample using the SN spectral fitting code, synapps, 



(Thomas, Nugent & Meza 2011). This code is based on synow 



Fisher 2000) and has the same assumptions, and so is limited to 



identification of features and estimates of their ejection velocities. 
It does not provide measurements of ion species abundances. The 
main advantage of synapps over synow is that it optimises automat- 
ically over the input parameters to perform a highly parameterised 
fit to the observed spectra. The outputs of synapps include ions con- 
tributing to the various spectral features and the velocities of the 
individual ions used in the fit. 



2.3 Optical photometry 

Multi-colour light curves were obtain ed at the robotic 2-m Liver- 
pool Telescope (LT; ISteele et al. 120041) located at the Roque de Los 
Muchachos Observatory on La Palma. The optical imager, RAT- 
CAM was used with gri fi lters, similar to t hose used in the Sloan 
Digital Sky Survey (SDSS; lYork et albOOOh . RATCAM has a field 
of view of 4.6' x 4.6' with a pixel scale of 0.135" (unbinned). 

Deep references images of the SN fields were obtained at 
>3 10 d post explosion at the LT for all the SNe listed in Table [2] 
on photometric nights. These were subtracted from the SN images 
at each epoch to remove host galaxy contamination. This subtrac- 
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Figure 2. Our HST NUV spectral sample. The top panel shows the SN la spectra entering the mean spectrum comparison, while the bottom panel shows the 
SNe that have been removed due to B - V colour cuts (PTFlOygu, SN 201 lek), along with those removed due to the effective phase cuts (-1.0 to +4.5 d). For 
each event, the effective phase is listed next to the SN name. The four SNe la in the sample for which B - V colours could not be calculated are excluded from 
the plot since colour corrections could not be applied. 
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Table 2. Log of photometric properties of the SN la sample. 



SN name 


MJD 


Effective 


B-band 


Stretch 


B-V 




B max 


phase (d)* 


max (mag) 






PTF09dlc 


55073.7±0.2 


2.8 


18.04±0.07 


1.05±0.03 


-0.16±0.05 


PTF09dnl 


55074.8±0.1 


1.3 


15.80±0.02 


1.05±0.02 


-0.02±0.01 


PTF09dnp 


55071.0 ±0.2 


5.8 


16.60±0.07 


0.96±0.05 


-0.1 5 ±0.06 


PTF09fox 


55131.7±0.2 


2.6 


18.25±0.06 


0.92±0.04 


0.00±0.04 


PTF09foz 


55131.8+0.3 


2.8 


17.81±0.10 


0.87±0.06 


0.03±0.08 


PTFlObjs 


55261.6±0.2 


1.9 


16.01±0.03 


1.08±0.02 


-0.09±0.02 


PTFlOfps 


55311.7+0.3 


10.6 


16.66±0.05 


0.73±0.03 


0.14±0.04 


PTFlOhdv 


55344.1 ±0.3 


3.3 


17.56±0.03 


1.05±0.07 


0.09±0.03 


PTFlOhmv 


55351.4+0.1 


2.5 


17.33±0.06 


1.15+0.01 


0.19±0.04 


PTFlOicb 


55360.6±0.1 


0.8 


14.48±0.03 


0.99±0.03 


0.06±0.02 


PTFlOmwb 


55390.7± 0.1 


-0.4 


16.81±0.04 


0.94±0.03 


0.03±0.03 


PTFlOndc 


55390.3±0.1 


5.8 


18.41±0.02 


1.05±0.03 


-0.05±0.02 


PTFlOnlg 


55391.5 ±0.2 


6.2 


18.60±0.03 


0.94±0.05 


0.13±0.03 


PTFlOpdf 


55407.4±0.3 


2.2 


- 


1.23±0.03 


- 


PTFlOqjl 


55418.9 ±0.1 


5.9 


17.79±0.03 


0.93±0.02 


-0.10±0.02 


PTFlOqjq 


55421.0±0.1 


3.5 


16.08±0.02 


0.96±0.02 


0.08±0.02 


PTFlOqyx 


55426.1 ±0.1 


6.8 


18.21±0.02 


0.85±0.01 


-0.10±0.02 


PTFlOtce 


55442.0±0.1 


3.5 


17.15±0.03 


1.07±0.02 


0.03±0.02 


PTFlOufj 


55456.5±0.2 


2.7 


18.31+0.09 


0.95±0.02 


-0.08±0.06 


PTFlOwnm 


55476.5±0.1 


4.1 


18.17±0.02 


1.01±0.03 


0.04±0.02 


PTFlOwof 


55474.2±0.1 


5.9 


17.93±0.07 


0.99±0.04 


0.10±0.04 


PTFlOxyt 


55490.9±0.2 


3.2 


18.40±0.04 


1.07±0.04 


0.19±0.03 


PTFlOygu 


55495.8±0.1 


-0.3 


17.29±0.04 


1.07±0.02 


0.44±0.03 


PTFlOyux 


55496 4+n 1 


7.1 


18.44+0.05 


0.83±0.01 


0.16±0.02 


PTFlOzdk" 












PTFlOacdh" 


55550.1 ±0.1 


9.1 




0.80±0.02 




PTFllkly 


55814.9±0.04 


2.8 


10.12+0.01 


1.00±0.01 


0.06±0.01 


SN20091e 


55164.9±0.5 


0.3 


15. 15 ±0.06 


1.08±0.01 


0.06±0.07 


SN2010ju 


55523.9±0.4 


5.5 


15.79±0.23 


1.05±0.58 


0.23±0.06 


SN2010kg 












SN2011by 


55690.6±0.1 


-0.3 


12.97±0.07 


0.93±0.02 


0.02±0.01 


SN2011ek 


55788.9±0.1 


3.7 


13.84+0.15 


0.90±0.02 


0.18±0.04 



"No multi-colour light curves were obtained for PTFlOpdf, PTFlOzdk, PTFlOacdh and SN2010kg. 
* Effective phase is the measured phase divided by the stretch. 



tion involves registering and combining five reference images with 
individual exposure times of 150 s to obtain a stack. For each SN 
image, the reference image was flux-scaled to it and then subtracted 
off using a PSF-matching routine. Where possible, the zero-points 
were calculated directly using SDSS stars in the fields of the SNe 
to calibrate to the SDSS photometric system. This was not possible 
for PTFlOnlg and SN20091e, which were not in the SDSS footprint. 
These were instead calibrated using standard star fields, obtained 
on photometric nights, observed before and after the deep SN ref- 
erence images. This is a preliminary calibration, and the final low-z 
sample of SNe that PTF will obtain will all be fully calibrated to 
SDSS Stripe 82 fields. 

B-band light curves for the SNe la sample were also obtained 
using the PTF search telescope, the P48. Pre-explosion reference 
images were available for the SNe discovered with PTF. The P48 
data were reduced by the Infrared Processing and Analysis Center 
(IPAC| pipeline (Laher et al. in prep.), and the light curves mea- 
sured using a custom-built pipeline, similar to the LT pipeline de- 
tailed above. The photometric calibration was performed to SDS S 
following a similar, but independent, method to lOfek et alj d2012h . 



Due to bad weather and scheduling difficulties, suitable 
LT+RATCAM light curves were not obtained for PTFlOpdf, 
PTFlOacdh, PTFlOzdk and SN 2010kg. For some SNe, it was pos- 
sible to supplement the light curves with data from the Faulkes 
Telescope North (FTN) on Haleakala, of the Las Cumbres Obser- 
vatory Global Telescop^l (LCOGT). These data have been reduced 
and calibrated in a similar manner to that described above for the 
LT data. Multi-colour light curves of PTF1 lkly (SN 201 lfe) were 
obtained with the Byrne Observatory at Sedgwick Reserve (BOS) 
0.8 m Telescope in g'r'i' band filters and calibrated to the SDSS 
system. 

2.4 Light curve fitting 

The optical gRri-bmd l ight curves were an alysed using the SiFTO 
light curve fitting code dConlev et al.l2008T) , which produces values 
for the stretch, maximum B-band magnitude, B-V colour at max- 
imum, and time of maximum light for each SN. These photometric 
parameters are listed in Table [2] SiFTO uses a time-series of spec- 
tral templates that are adjusted to recreate the observed colours of 



2 http://www.ipac.caltech.edu/ 3 http://lcogt.net/ 
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Figure 3. The mean spectra of the HST sample are shown for blue and red SNe la as the solid lines in these plots. The left panel shows the SNe before 
colour-correction, when split based on the B — V colour obtained from light curve fits, with blue (B -V < 0.025) and red (B — V > 0.025). In the right panel, 
the CCM-colour corrected SN mean spectra are shown. The mean spectra were constructed as described in Section |3. 1| The shaded regions represent the 90 
per cent confidence levels from a bootstrap resampling. 



the SN photometry at each epoch, while also adjusting for Galac- 
tic extinction and redshift (i.e., the ^-correction). The P48 B-band 
data is not used in the B - V colour measurements, just for con- 
straining the date of maximum and the stretch. Provided there are 
measured colours for a SN, SiFTO can also interpolate to obtain the 
peak magnitude in a chosen rest-frame filter. The magnitudes and 
colours quoted throughout the paper have been measured at B-band 
maximum. The colour output of SiFTO is C, a weighted average of 
the U — B and B—V colours at B-band maximum. However, because 
our bluest filter is the g-band, SiFTO C and B—V are equivalent 
for our dataset. Therefore, we will use the term B—V throughout 
to mean the C colour a t B-band maximum light. Tests have been 
performed in Guv et al. (2010) comparing SiFTO and the popular 
SALT2 l lGuv et alj|2007l) light curve fitter, and good agreement is 
found between their outputs, despite the differences in their meth- 
ods. The effective phase values quoted in Table[2]are the measured 
phases divided by the stretch as measured from the light curve fits. 
A comparison analysis has been performed using the phases before 
correcting for stretch and does not change the results. 



2.5 Colour correction of UV spectra 

Dust extinction towards SNe la can affect the measured 'colour', 
with different amounts of extinction resulting in varying colours 
being observed. If the intrinsic colours of SNe la were assumed to 
be the same for all objects and Milky Way-like dust was responsi- 
ble for the observed variations in B - V colour, then when using the 



Milky Way dust extinction CCM law dCardelli. Clayton &~M athis 
1989), a wavelength parameterisation with a selective-to- total ex- 
tinction value of R B = 4.1 (R v = 3.1) should be applicable. How- 
ever, it has been found that when the value of R B i s allowed to vary, 
an R B a 3 is favoured by SNe la dTriprjll998l: iTripp & Branch! 



ll999tlConlevetalj2007h . This suggests that either the dust in SNe 
la host galaxies is not consistent with Milky Way-like dust, or that 
there are intrinsic variations in the colours of SNe la, similar to the 
observed variation in light curve stretch. Here we correct for the ob- 
served B—V colour variations, to remove the differences in B - V 
colour (dust and intrinsic colour) that could affect the measured flux 
in the spectra. 

We investigate d two ways to ap ply the colour correction, us- 
ing either the SALT dGuv et al.l2007h or the Milky Way dust CCM 
colour law. The SALT and CCM colour laws have similar forms at 
optical wavelengths, but are different in the NUV, where the SALT 
(SALT1, SALT2) colour laws have steeper slopes. Both the SALT 
colour laws were applied using the SiFTO C, corrected to a SALT 
colour, using the conversion relation from iGuv et al.l d2010l) . The 
chosen value of total-to-selective extinction does not matter since 
the spectra are normalised in the region 4000-4500 A and measure- 
ments of wavelength positions are also not affected by this value. 

Fig. [3] shows the low-z sample split into two bins at the po- 
sition of the mean C of the sample, before (left panel) and after 
(right panel) a CCM colour-correction was applied. Before correc- 
tion, SNe with bluer C have bluer spectra, as is expected. After 
CCM-correction, there is better flux agreement between the C split 
samples. The CCM colour law correction gave a more consistent 
agreement between spectra than either the SALT1 or SALT2 colour 
laws, which overcorrected the colour of the NUV spectra making 
the redder SNe significantly bluer than the originally bluer SNe. 
E08 used the SALT1 law because of better post-correction agree- 
ment between the low and high C SNe, but we note that in fig. 8 of 
E08 the difference between the CCM and SALT1 colour law cor- 
rections is small. Therefore, for consistency the CCM colour law 
has been used to correct both the low-z and intermediate-z samples. 
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Figure 4. Comparison of the mean spectra of the low-z sample (17 SNe la) to that of the intermediate-z sample (15 SNe la) in the top panel. The shaded 
regions represent the 90 per cent confidence levels from a bootstrap resampling. The region through which the normalisation was applied is marked, along 
with positions of the A[ and A2 NUV features. The lower panel shows the percentage dispersion at each wavelength (bootstrap resampling error as a percentage 
of the flux of the mean) for the HST and Keck samples respectively. An increase in the dispersion in the NUV compared to optical wavelengths is seen. A 
peaked appearance is seen in the dispersion of both samples around ~3700 A, which is found to be caused by an increased dispersion in the blue wing of the 
Can H&K feature. 



2.6 Host galaxies 

The redshift of the host galaxies were obtained in order of prefer- 
ence i) from host galaxy lines identified in the optical SN spectra, ii) 
from galaxy spectra obtained from the NASA/TPAC Extragalactic 
Database^ (NED) or iii) from the SDSS (labelled 'NED' or 'SDSS' 
in Table [TJ. However, it was not possible to obtain a host galaxy 
redshift for two SNe, PTFlOufj and PTFlOqyx (labelled 'template' 
fit in Table [T}. A faint host at the position of PTFlOufj was de- 
tected in a 200 s ^-band image obtained with Keck+LRIS, but a 
1800 s combined spectrum, obtained with KECK+LRIS using the 
600/4000 blue grism and 400/8500 red grating, did not reveal any 
identifiable features that could be attributed to a host galaxy. For 
PTFlOqyx, no host was detected in a 200 s ^-band image obtained 
at KECK+LRIS to a limiting magnitude of ~25.5 mag. 

For these two SNe, redshift s were obtaine d using the template 
matching routine superfit (Howell et a I]|2005h by selecting the av- 
erage redshift obtained from all the available spectra for these SNe. 
This redshift was then used to correct the spectra to the rest -frame. 
These SNe were excluded from the wavelength shift and velocity 
analysis in Section [4] because of the larger uncertainties in their 
redshifts. 

The host galaxy M slc i lar of our SN la sample were calculated 
4 http://nedwww.ipac.caltech.edu/ 



using the PEGASE code fce Borgne & Rocca- Volmerangeir2002l). 
with S DSS multi-band photometry and the method of Sulliv anet al.l 
(2010) to fit the SED of the host photometry to a series of galaxy 
templates. The values of these fits are listed in Table[T]and the dis- 
tribution of the host MjteUaj is shown in Fig.Q] For SNe in the sam- 
ple that do not have SDSS photometry available, the magnitude of 
the galaxy was estimated using the gn'-band deep reference images 
obtained with LT, from which the M ste iiai could be estimated. 

If the host galaxy was not detected (PTFlOqyx, PTFlOufj), a 
M 6te uar of 10 6 M Q was assigned. However, there are three galaxies 
at separations of 1 - 5' from the position of PTFlOqyx, that may be 
potential host galaxies. The nearest potential host galaxy, VVDS 
020218384 is at a separation of 0.9"and has a redshift of 0.068, 
within the uncertainti es of the SN redshift value of 0.065+0.005 
dLe Fevre et all [2005 ) . This suggests that the host of PTFlOqyx 
may not be very faint/low mass but more massive and distant. How- 
ever, PTFlOqyx are not included in the wavelength shift and veloc- 
ity studies, and so will not affect our conclusions. 

2.7 Hubble residuals 

For cosmology using SNe la, the peak fi-band magnitude is 'cor- 
rected' using empirical relation between light curve stretch, light 
curve colour, and host M slc ii ai .. The residuals on the Hubble diagram 
after correction for some or all of these parameters are studied to 
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Figure 5. The SN la low-z subsample split based on stretch, with a cutoff between high (blue) and low (red) stretch of 1.01. The bottom panel shows the 
dispersion of the high- and low-stretch means. The numbers in parentheses are the number of SN la spectra that went into each mean spectrum. 

(e.g. ISullivan et alj|200l ICooke et al.ll20lT| ; lFolevll2012h . and to 
understand whether the UV dispersion originally discussed by E08 
is an intrinsic property of SNe la or connected to this claimed evo- 
lution. The impact of these trends on the cosmological utility of 
SNe la can be gauged by examining possible correlations between 
the NUV spectra and the photometric residuals on the Hubble di- 
agram. We compare low and intermediate-z samples, along with 
comparisons of the HST SNe split based on stretch and Hubble 
residuals. The spectra split by B - V colour have been described in 
Section |2~5l when discussing the colour corrections, and so are not 
described again here. 



3.1 Constructing mean spectra 

To look for any potential evolution with redshift in the SNe la sam- 
ple, the new low-z sample must be compared to a sample of SNe la 
NUV spectra at higher z for which we use the intermediate-z sam- 
ple of E08. The spectra from this intermediate-z sample have been 
analysed in the same manner as our low-z sample, and we have 
colour-corrected them using the CCM colour law, as described in 
Section l231 For consistency, we have also refitted their light curves 
using the same version of Si FTO that is used for our low-z sample. 
Gonzalez-Gaitan ~et alj d201lf) investigated the B - V colour-stretch 
relation of the SNLS SN la sample and found no correlation be- 
tween light curve stretch and B - V colour for 'normal' SNe la. We 
also find no correlation between stretch and B - V colour for our 
low-z sample. 

The mean spectra are created by rebinning the input spectra to 



look for possible correlations with spectral features. These could be 
used to further reduce the scatter in the Hubble diagram or provide 
alternative correction methods. 

The residuals, m^ sW are calculated using 

mT d = <" -M B -5 lo gl0 D L -25+ a(s - 1) - BC (1) 

where m™ v is the measured B-band magnitude, D L is the luminos- 
ity distance in units of Mpc, s is the stretch and C is the colour of 
the SN. a, B and M B are the coefficients of the best- fit cosmologi- 
cal par ameters to the data following the technique in lSullivan et al.l 
J201 1). We assume a standard flat-ACDM cosmology (fl„, = 0.27, 
fl A =0.73, w = -1). Residuals can also be calculated removing the 
correction for stretch, 

m r^.no , = mT _ Me _ 5 hgw Dl _ 25 _ pC (2) 

These residuals can be used to look for correlations with the wave- 
length/velocities of NUV features. We apply a redshift cut of z > 
0.02 (to exclude SNe that are affected by host galaxy peculiar ve- 
locity components), which excludes seven SNe, leaving 1 1 SNe for 
which the Hubble residuals can be calculated (after the phase and 
B —V colour cuts detailed in Section[3~lt. 



3 RESULTS I - MEAN COMPARISONS 

This section details comparisons performed between the mean 
spectra described in Section |3~T1 We seek to test the claim that the 
NUV spectra evolve with redshift as suggested by previous works, 
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3 A per pixel, and then normalising the spectra through a box fil- 
ter over a wavelength range 400CM-500 A. A three-sigma clipped 
mean is used and the error on the mean is determined by bootstrap 
resampling. In the mean plots shown here, the dashed lines refer to 
lower and upper 90 per cent confidence limits. The results are not 
changed by varying the normalisation region in the optical wave- 
length region. 

SN la spectra in a suitable phase range must be selected when 
constructing the NUV mean spectra so that phase-related effects do 
not unnecessarily influence the mean spectrum comparisons. We 
have chosen an effective phase range of -1.0 to +4.5 d for both the 
low-z sample and the intermediate-z sample from E08. Coo ke et"aL1 
j201 ll) chose an effective phase range of -0.32 to +4 d, similar to 
our chosen values. We have also placed cuts on the stretch (0.7 < 
s < 1.3), and B - V colour (-0.25 < c < 0.25) of the SNe in both 
samples. T hese are the same c uts as applied in the cosmological 
analysis of lConlev et al.1 J2011I) . No SNe la in our low-z sample are 
removed by either of the stretch cuts or the lower B - V colour or 
phase cut. We lose nine SNe from the low-z sample due to the upper 
phase cut and one SN (PTFlOygu) because of the upper B- V colour 
cut. SN 201 lek is also excluded, since it displays an unusual NUV- 
B colour of 1.0 mag, which will be discussed in Section l4~4l This 
leaves a total of 17 SNe to be used in the mean low-z spe ctrum. The 
preliminary low-z mean spectrum of ICooke et al.1 d201 lb contained 
10 HST SNe la and one historical SN la from the literature. From 
the initial intermediate-z sample of 33, 17 SNe are removed due to 
the phase cuts and one SN due to the B - V colour cuts, leaving 15 
SNe for the intermediate-z mean. 

Choosing SN la spectra in a suitably small effective phase 
range, yet that still contains a statistically significant number of ob- 
jects, will limit the effect of phase variations when making compar- 
isons between mean spectra. Despite this, phase variations within 
this range will still occur; for example, SNe cool with time, and 
their line velocities and NUV flux decrease. However, the effect 
of phase variations can be minimised by ensuring that when com- 
paring mean spectra, the phase distributions (along with the stretch, 
B—V colour and host galaxy distributions) of the samples are drawn 
from the same parent populations. As described in Section |2~T1 K-S 
tests were performed and for all these properties, there is no indi- 
cation they are drawn from different parent populations. 

3.2 Evolution in the mean SN la spectrum with redshift 

The mean spectrum at low-z is compared to that of the 
intermediate-z Keck SN la sample from E08 in Fig. [4] This is the 
first time that a large sample of low-z SN la NUV spectra with 
colour corrections applied has been available for comparison with 
higher-z samples. As described in Section |3~T| the effective phase, 
stretch, B - V colour and host M stc n ar distributions of the two sam- 
ples are found to be well-matched. The mean spectra of the two 
samples are seen mainly to agree with the uncertainties (90 per 
cent confidence limits). However, at shorter wavelengths (< 3300 
A), the low-z mean spectrum has less flux than the intermediate-z 
sample. 

To investigate this further, a random number of spectra (8- 
15) were used to make the low- and intermediate-z mean spectra 
and this was repeated 20000 times. We compute the flux through 
a 'UV box filter, in the region 2900-3300 A, for each spectrum, 
and find that 99.8 per cent of the time (3.1-cr), the intermediate-z 
mean spectrum ha s a greater flux than the low-z mean spectrum. 
ICooke et al.l J20T l|) also found a lower flux at shorter NUV wave- 
lengths when compared to the sample of E08. The choice of this 



wavelength range o f 2900-3300 A was motivated by the results of 
IWalker et al.1 J2012h . which showed using radiative transfer mod- 
elling that this is the region where the most variation due to metal- 
licity effects is expected. The importance of this flux evolution with 
redshift will be discussed in Section l4"l4l 

We also compare in Fig.|4]the dispersion in the mean low- and 
intermediate-z SN samples. The dispersion is calculated as the flux 
error as a percentage of flux. There is an increased dispersion at 
wavelengths < 3700 A, as found in previous studies. A peaked fea- 
ture is visible at ~3500-3800 A in the dispersion of both samples. 
This is found to be caused by an increase in the bootstrap resam- 
pling error in the lower part of the blue wing and the minimum of 
the Ca 11 H&K feature and the relatively small flux in this region. 

We investigated the effect of host galaxy contamination on 
the results obtained for the mean spectrum of both the low- and 
intermediate-z sample. As described in Section l2?2l the HST sam- 
ple is free from host galaxy contamination due to the narrow slits 
used to obtain the spectra. The host galaxy contamination of the 
intermediate-z sample was estimated in E08 and the continuum 
flux levels in the NUV were found to be well constrained using 
the multi-band u'g'r'i'z' SNLS photometry. The removal of host 
galaxy emission lines from the spectra of this sample was more 
complicated and is the cause of the discrepancy in the trough of the 
Can H&K feature, but this does not affect the overall flux levels. 
The effect of the colour correction used on mean spectra is also in- 
vestigated, and the flux offset between the low- and intermediate-z 
samples remains when different colour laws are used. 

3.3 Mean spectra split by light curve shape 

Our sample can be split in various different ways to investigate 
possible astrophysical effects in SNe la. One of the most obvious 
is to split the sample according to the SN light curve stretch. We 
have split the sample into bins at the position of the mean stretch 
(s = 1.01) of the sample. This results in 9 and 8 SNe in the low- 
and high-stretch bins respectively. Fig.|5]shows a comparison of the 
low- and high-stretch mean spectra for the low-z sample. A blue- 
shift in the wavelength positions of the NUV features is clearly 
seen for the higher stretch SNe. In particular, the blue wings of Ca 11 
H&K feature at -3700 A and the NUV features, lj and A 2 at 2920 
and 3180 A respectively, are blue-shifted in high-stretch compared 
to low-stretch SNe la. We have shown that the phase distributions 
of the low- and high- stretch samples are well-matched (97 per cent 
probability they are from the same parent populations). Therefore, 
the effect of different phase distributions can be ruled out as the 
cause of the velocity shift seen in Fig. [5] 

Fig. [5] shows that the low-stretch SNe in our sample have 
slightly more NUV flux at the positions of the A\ and A 2 features. 
As will be discussed in Section [5] this is most likely linked to the 
velocity differences seen in the NUV, with higher velocities causing 
more line blanketing and hence, a greater suppression of the NUV 
flux for high, compared to low velocity events. 

The dispersion of the two stretch samples is also shown in Fig- 
ure|5] where an increased dispersion towards shorter wavelengths is 
observed. This is similar to the increased dispersion found in both 
the low and intermediate-z mean spectra. Similarly to the spectra 
split by redshift, a peaked feature is present at a slightly shorter 
wavelength than ~ 3700 A, when the sample is split by stretch. The 
peak of the high-stretch dispersion is shifted by 21 A with respect to 
the low-stretch dispersion. This feature is caused by a combination 
of an increased bootstrap resampling error at the minima and lower 
half of the blue-wings of the Can H&K features in both spectra, 
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Figure 6. The mean spectra of HST SN la subsamples split by uncorrected Hubble residual value (left) and by corrected Hubble residual value (right). The 
division value is a Hubble residual of zero. The numbers in parentheses are the number of SN la spectra that went into each mean spectrum. 



and the smaller flux in this region than at the blue-edge continuum 
of the feature. 

3.4 Mean spectra split by Hubble residuals 

One interesting parameter upon which our sample can be split is 
Hubble residual values. These values were calculated as described 
in Section |2~7l In Fig.|6]the mean spectra split by uncorrected Hub- 
ble residual and by corrected Hubble residual are shown. The un- 
corrected Hubble residuals are the values before corrections for ei- 
ther light curve stretch or colour are applied (see Section |2~7l l. The 
corrected Hubble residuals are those after these corrections have 
been applied. The velocity shift in the blue wing of the Can H&K 
feature for SNe la with negative Hubble residuals seen in the left 
panel of Fig.|6]is in agreement with the results of the mean spectra 
split by stretch. There is good agreement between the two spectral 
comparisons shown in Fig. [6] with no obvious differences between 
the positive and negative Hubble residual samples. 



4 RESULTS II - ANALYSIS OF INDIVIDUAL SPECTRA 

Although the comparison of mean spectra is a useful tool for assess- 
ing evolution in SN properties as a whole, the analysis of individ- 
ual spectra can be used to investigate colour variations, along with 
correlations between spectral properties and other SN la proper- 
ties. Individual line measurements are also advantageous over mean 
spectrum comparisons because it is possible to correct for phase 
variations, which can occur even within the narrow phase window 
used for the mean comparisons (-1.0-4.5 d). The same cuts on the 
sample as in Section[3]are used. 

The NUV colours can be calculated from the spectra and com- 
pared with the B - V colour obtained from the light curve fits to 
look for increased dispersion at NUV wavelengths. An analysis of 



the features in individual spectra is also performed, including the 
wavelength shifts of the Si n 4130 A, Can H&K, A 2 and Ai features. 
The At and A2 features are defined as the NUV peaks at ~2920 and 
~3 180 A respectively as shown in Figure|4] These features are sug- 
gested from modelling to be depend ent on the metallicity of the 
progenitor and its com position (e.g. iLentz et alJlioOOl : ISauer et al.l 
2008]; IWalker et al.ll2012h . Therefore, a shift in these features may 
be driven by the properties of the progenitor system. 

In Fig. [7] the individual spectra that are used in the measure- 
ments presented in this section, are shown for the Ca 11 H&K and 
Sin 4130 A features on the left, and Ai, A 2 features on the right. 
The SNe spectra are further split based on their stretch, into low- 
and high-stretch bins at a position of s = 1.01. It can be seen clearly 
that there is a stretch-dependent shift in the wavelength position of 
the Ca n H&K, A 2 and Si 11 4130 A features, with higher stretch SNe 
having greater velocities, similar to the results of Section |33l These 
trends between the wavelength shift of the features and stretch will 
be investigated further in the rest of this section. The increased 
dispersion at shorter wavelengths is also clearly visible in Fig. [7] 
where there is a much larger variation in the positions of the A 2 , 
and particularly Ai features compared to the relative homogeneity 
(apart from velocity differences) of the Can H&K feature. This in- 
creased dispersion will investigated further using synthetic filters 
in Sectionl44l 

4.1 Line measurements 

To measure the position of the minima and maxima of the vari- 
ous spectral features, we fit a Gaussian to the line profile. This is 
done by subtracting off the pseudo-continuum, then choosing max- 
imum values o n each side of the feature and performing a series of 
IDL's MPFIT iMarkwardtll2008h fits around this initial value. The 
positions of the pseudo-continuum points are then varied within a 
range of +30 A and the minimum wavelength is chosen as the min- 
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Figure 7. The individual SN la spectra that are used in producing the mean spectra presented in Section [5] are shown for the Can H&K and Si n 4130 A 
features on the left and Xi, Aj features on the right. The SNe spectra are also split into high-stretch (blue) and low-stretch (red) bins. The numbers in the 
legends represent the number of SN spectra plotted for each stretch bin. The spectra have been smoothed with a boxcar average of 1 1 pixels. The spectra have 
been normalised in the region 3400-4300 A for the left panel and 3000-3400 A for the right panel. 




Figure 8. The measured velocities (corrected for phase variations) of Can H&K and Sin 4130, along with the wavelength of the Ai feature are plotted as 
functions of each other. The left panel is Ai against Can H&K, the middle panel is Can H&K plotted against Sin 4130 A and the right panel is Az against Sin 
4130 A. The high-stretch (s > 1.01) and low-stretch (.v < 1.01) samples have been displayed in blue filled circles and red open circles respectively. There is 
one less SN in the high-stretch bin for the plots containing Ai, because the A2 feature could not be identified for one SN (PTF09dlc). The best fit lines to the 
data in each panel are plotted as solid, black lines. The significances of a non-zero gradient are 3.4, 2.8 and 2.2-cr from left to right. 



imum of the fit with the lowest x 2 value. The wavelength position 
errors are estimated from the range in these fit values. The total un- 
certainties of the wavelength measurements include the uncertainty 
in measuring the position, selecting the pseudo-continuum and the 
redshift uncertainties. All the spectra were visually inspected to en- 
sure the fit and the chosen pseudo-continuum values are sensible. 

iFolevI J2012h (whose data will be compared to ours in Sec- 
tion !4.5t measure t he line position of the Ca 11 H&K feature follow - 
ing the method of iBlondin et alj 120061) and the culling criteria of 
iFolev. Sanders & Kirshner] J201 lb - This method involves smooth- 



ing the spectra using an inverse-variance-weighted Gaussian filter. 
The smoothed spectra are resampled onto a fine wavelength scale 
of 0. 1 A and the minima of the features are measured. If two min- 
ima are identified in the Ca n H&K feature, they are classified sepa- 
rately as "blue" or "red" compon ents. Five SN spectra present in the 
Keck sample are also part of the lFolevI d2012l) sample. We compare 
o ur measureme nts of the Can H&K velocities to those presented 
in IFolevI J2012h and fi nd our measured values are lower than the 
values of Folev ( 2012) by values ranging from ~650-3800 km s -1 



iBronder et alj d200 8) also used a Gaussian fit to the entire feature, 
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while ISilverman et al.l J2012bh fitted the entire Ca n H&K line pro- 
file using a cubic spline to identify the position of the minimum. 

We have chosen not to follow the method of Foley (2012) 
since the presence of additional features present in the Ca n H&K 
features could bias the measurements of these parameters. There- 
fore, we fit a Gaussian to the overall Can H&K feature. Although 
this will include a contribution from the Sin feature, which has 
been identified in the Can H&K feature using the spectral fitting 
code, synapps, the contribution from Sin can be estimated using 
the velocity outputs from synapps (see Section|4~i 



4.2 Correction for phase variation 

As for our mean spectrum comparisons in Section [3] we measure 
wavelength positions/velocities for the studied features (Ca n H&K, 
A.I, A 2 and Sin 4130 A) in the SN spectra in the effective phase 
range of -1.0 to +4.5 d. However, variations as a function of time 
within this chosen phase range may still be present and must be 
corrected for. We correct these variations by fitting the wavelength 
positions/velocities as a function of phase for our SN samples and 
adjusting to a phase of d. We simultaneously fit the values of both 
the HST and Keck spectra for the Can H&K velocity, X% wave- 
length and Sin 4130 A velocity, while for the A { feature, we fit just 
the HST spectra due to their higher signal-to-noise (S/N). The fits 
are performed using equations (depending on whether wavelength 
or velocity is being corrected) of the form (here for velocity), 



I'D 



(3) 



where v, is the velocity (kms~') measured from the spectra, Vo is 
the velocity (km s -1 ) at maximum, t is the effective phase (d) of the 
spectrum and v is gradient of the fit ( kms^'d -1 ). The values of the 
gradient of the fits are shown in Table [3] along with the maximum 
correction to d that was applied. The NUV-optical colours mea- 
sured from the NUV spectra are also corrected for phase variations 
following the same method. For the NUV-optical colour measure- 
ments, we fit the values of only the HST sample. The measured (and 
phase-corrected) wavelengths and velocities are given in Table lAll 
Ideally the velocity evolution of the individual SNe la should 
be studied but given the single epoch spectra available for our 
sample, this was not attempted. However, a decrease in veloc- 
ity/increase in wavelength is seen as a function of time for our 
sample, as has been found using meas urements of ind i vidual 
SN spectra as a function of t i me (e.g. [Hachinger et al.l 120061 : 
iFolev, Sanders & Kirshnerll201ll : ISilverman et al.l l2012bl) . There- 
fore, on average, our corrections should remove this phase depen- 
dence in the range of -1.0 to +4.5 d. Our NUV HST Cycle 18 
data, with multi-epoch spectra of 4 'normar SNe la will provide 
improved measurements of the diversity of the phase evolution of 
NUV features for future studies. 



4.3 Ca 11 H&K, Si 11 4130 A and A 2 positions 

Fig.[8]shows the wavelength positions of the Can H&K, Sin 4130 
and A 2 features plotted as functions of each other. The left panel is 
A2 against Can H&K, the middle plot is Can H&K plotted against 
Sin 4130 A, and the left plot panel is A 2 against Sin 4130 A. The 
high- and low-stretch SN samples are also shown, split at a value of 
s = 1.01. The best fit lines to the data for each comparison are also 
measured and have significances from zero of 3.4, 2.8 and 2.2-cr 
for A 2 against Can H&K, Can H&K plotted against Sin 4130 A 
and A 2 against Sin 4130 A, respectively. 



Table 3. The gradient of the phase corrections that were used to correct the 
wavelengths of the studied feature and colour to a phase of B-band maxi- 
mum, and size of the correction applied. 



Parameter 


Gradient 


Size of 






correction 


Can H&K 


-280±230kms-'d- 1 


< 1280 km s -1 


h. 


6.12±2.61 Ar 1 


< 25 A 


M 


-1.64±2.75 Ad~' 


< 7 A 


Sin 4130 


-a^Okms-'d- 1 


< 372 km s" 1 


UV-b 


0.06±0.02 mag d" 1 


< 0.3 mag 
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Figure 9. The UV-b colour (after phase and SiFTO B-V colour correction) 
is shown as a function of B-V colour. The outlier in the plot is SN 201 lek, 
which has a very red UV-b, along with one of the reddest B-V colour 
values. 



The Sin 4130 A absorption feature has been previously con- 
sidered a good measure of the photospheric velocity and there is a 
correlation between this feature and the Ca 11 H&K and A 2 features. 
However, the Sin 4130 A feature is investigated using synapps and 
it is found that S 11 and Co n are found to contribute in this region. 
This could explain the lack of a stronger correlation between the 
measurement of the Sin 4130 A 'notch' and the Can H&K ve- 
locity seen in Fig. [8] since this is not measuring Si n alone. Fur- 
ther investigation of the NUV features using both correlations with 
photometric properties, and using synapps spectral fitting will be 
presented in the following sections. 



4.4 NUV-optical flux comparison 

The evolution with redshift in the NUV spectral features, as well 
as the increased dispersion at NUV wavelengths shown in Section 
|3.2| can also be investigated using synthetic filters. As discussed in 
Section [3~2l to investigate the significance of the observed evolution 
with redshift, we define a 'UV box filter between 2900-3300 A, 
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since this is the spectral region which is shown by Walk er et al.l 
d2012h to have the largest metallicity variations. This 'UV' filter is 
100 A redward of the 'UVT filter used in E08. We compare the 
flux through this filter to that in a 'b'-band filter with a wavelength 
range of 4000-4800 A. 

We investigated the size of the evolution with redshift in the 
mean spectra by comparing the mean weighted UV - b colours of 
the low- and intermediate-z samples, and find that they have values 
of -0.18±0.03 and -0.36±0.04 mag, respectively. This difference 
of 0.18 mag (3.6-cr effect) is in agreement with the results of the 
mean spectrum comparison at different redshifts, shown in Section 
13.21 where the low-z is redder at NUV wavelengths compared to 
the intermediate-z sample. We discuss the significance of this in 
Section|5] 

Fig. [9] shows the measured UV - b colours as a function of 
B-V colour (after applying the CCM-law colour correction and the 
phase correction given in Table[3}. The same phase and B-V colour 
cuts as used for creating the mean spectra are used for this sample. 
SN 201 lek was excluded from the mean spectrum analysis because 
of its unusually red UV - b colour that is a clear outlier (UV - b of 
1.0 mag) when compared to the SNe shown in Fig. [9] We find no 
correlation between the UV - b colour and stretch. However, the 
scatter in U V - b is significantly larger than in B-V, with an r.m.s. 
of 0.21 versus 0.08 mag (excluding SN 201 lek). This is consistent 
with previous studies that noted a larg er dispersion in N UV spectra 
compared to the optical (e.g. E08; iFolev et alj|2008al) . Again, we 
discuss the possible origin of these effects in Section|5] 

4.5 Ca ii H&K velocity 

The Can H&K wavelength positions for both the HST and Keck 
sample were measured and corrected for phase variations as de- 
tailed in Sections |4. 1 1 and |4~2l These wavelengths were then con- 
verted to velocities using the relativistic Doppler formula and a 
rest wavelength taken as 3945 A (the gf-weighted average of the 
3933, 3969 A features). The velocities of the Can H&K feature 
for the HST and Keck samples are shown in Fig. [10] as a function 
of stretch in the left panel and B-V colour in the right panel. 
lFolevN2012l) have also measured the Ca n H&K velocity of a sam- 
ple of SDSS and SNLS SNe la. Their measured values which are 
within our chosen phase range are corrected to maximum light us- 
ing our phase relationship and are also plotted in Fig. [10] No corre- 
lation between stretch and B-V colour is app arent in our sample , 
in agreement with previous larger studies (e.g. iBrandt et alj|20ld : 
iGonzalez-Gaitan et alj|201 IT) . 

A trend of increasing velocity with increasing stretch is seen in 
the HST and Keck samples. A linear fit was performed to the HST 
sample and is shown in Fig. [10] The fit has a significance from zero 
of 3.4-cr. This higher velocity with stretch trend seen for the HST 
sample is consistent with Fig. [5] where the Ca n H&K feature is 
shifted to the blue for the higher stretch SNe. A contribution to the 
Can H&K region from Sin is identified using the spectral fitting 
code, synapps. The Si n component is weaker in all cases than the 
contribution from Can. The identified trend in the observed data 
between Can H&K velocity and stretch can be investigated using 
the velocity outputs from synapps. A correlation between the fitted 
Can velocity and stretch is found at a significance of 2.4-cr from 
zero. The velocity of the Si n component is also found to correlate 
with stretch (2.3-cr). This suggests that the identified trend of in- 
creasing velocity of the overall Ca n H&K features with increasing 
stretch (luminosity), is a result of a trend in both in the Ca n and 
Si n components. 



In the right panel of Fig. [10] the Ca n H&K velocities as a 
function of B - V colour are shown, for the HST, Keck and Foley 
samples. No trend with B-V colour is id entified in eit her the HST 
or Keck samples, while for the sample of Foley (2012), their iden- 
tified trend of increasing col our with increasing Ca n H&K velocity 
is clearly seen in their data. Silve rman et al.l d2012ah also found no 
trend between Ca n H&K velocity and B-V colour for their low-z 
sample. 

The sample of IFolev! d2012l) includes only SNe la with s< 
1.05, while the HST data contains 7 SNe with stretch values greater 
than this. When these high-stretch SNe are excluded from our low-z 
sample, the trend between Ca n H&K and stretch drops to less than 
l-<x from zero. Analysis of the full range of stretch values used in 
current cosmological samples demonstrates the trend of increasing 
Can H&K velocity with stretch. The reason for the observed dis- 
crepancies between Ca n H&K velocit y and B-V colour, when 
comparing the HST sample, to that of Foley! d2012h is unclear. A 
possible reason may be the differences in the velocity measure- 
ment techniques used, as discussed in Section l4~T1 or the definition 
of B - V colour. 

The host galaxy properties of SNe la have been demonstrated 
to have a significant impact on the properties of the SNe they host. 
We investigate here if there is a relation between host galaxy M stc n al 
(a good proxy for metallicity) and Ca n H&K velocity. We find a 
correlation between these two parameters for the HST sample at a 
signific ance of 1.7-cr from zero. This is in agreement with the re- 
sults of lFolevi d2012l) , who also find a similar correlation in their 
sample. The left panel of Fig. [TT] demonstrates this relation, show- 
ing th e velo c ity of the Can H&K line of the HST and the sam- 
ple of iFolevI d2012h as a function of the log(M stc n al ). This correla- 
tion also agrees with the relation between Can H&K and stretch 
seen in Fig. [10] since it has been shown that less massive, star- 
forming gala xies preferentially host SNe la with broader, b righter 
light curves dHamuv et al1 [l996. 2000; lSullivan et alj|2006l) . From 
Fig- El it can also be seen that there is a lack of low velocity (low 
stretch) SNe in hosts with low M stc n ar . 

To investigate if the Can H&K velocity and Miliar correlation 
is caused by the previously identified stretch and M^na,. correlation, 
we remove the observed Ca n H&K velocity and stretch correlation 
shown. In the right panel of Fig.QT] the Can H&K velocity (after 
correlation for the Can H&K velocity versus stretch has been re- 
moved) as a function of M stc n ar is shown, and the relation between 
Ca ii H&K velocity and Miliar is found to be removed. This shows 
that the observed Can H&K velocity versus M sli: i lar correlation is 
driven by the well-known correlation between host galaxy M slc i lar 
and light curve stretch. 

4.6 Wavelength of NUV A 2 feature 

The emission peak corresponding to the A 2 was chosen by selecting 
the closest feature to the A 2 of 3180 A as defined in E08 and is 
shown in Fig [4] Fig. 1 121 shows the position of the A 2 feature as a 
function of stretch. A clear correlation between the position of the 
A 2 feature and stretch is seen: lower stretch SNe la have a redder 
wavelength position of the A 2 feature than higher stretch SNe. A 
best linear fit is determined with a significance of 6.2-cr from zero 
for the combined HST and Keck samples. Therefore, this trend is 
even stronger than that identified for the Can H&K velocity. We 
find no trend between the A 2 wavelength and B-V colour. 

The position of the A 2 feature as a function of M slc u ar shows 
a trend of increasing blueshift with decreasing host M stc u ar , with a 
significance of 2.0-cr. More generally, there is a lack of SNe la with 
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Figure 10. The wavelength position of the Ca n H&K feature for the HST (blue, solid circles) and Keck (red open circles) samples against stretch is shown in 
the left panel, while the right panel shows the velocity against B — V col our. The Ca n H&K velocities have been phase-corrected using the relation in Section 
|4.2| The grey, open diamonds are the measured Ca n H&K values from Foley (2012), which we have corrected using the same relationship as for the HST and 
Keck data. The solid blue lin e is the best fit to the HST data with a significance of 3.4-cr from zero. The same sample cuts as for the HST and Keck samples 
have been also applied to the Foley (2 0121) data. 
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Figure 11. The measured velocities (corrected for phase only) of the Ca n H&K feature for the HST sample and the sample from Foley (2012) are shown as 
a function of log(M ste n ar ) in the left panel. The solid blue line is a linear fit to the HST data only, showing a correlation between Can H&K velocity and host 
galaxy M stc iiar (17-<x from zero). iFolevI | |2012|) also identified a similar trend in their data. The right panel shows the Ca n H&K velocities for the samples 
after the correlation between Ca 11 H&K velocity and light curve stretch has been removed. A linear fit to the HST data is plotted as a solid blue line and the 
correlation between Ca 11 H&K velocity and host galaxy Miliar is seen to be removed. 
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Figure 12. The wavelength position of the A 2 feature for the HST and Keck 
samples, as a function of stretch is shown. A best linear fit to both data sets 
is plotted as a solid, black line in the left panel and has a significance of 
6.2-cr. A linear fit to the HST sample alone has a significance of 3.7-cr. 



large A 2 blue-shifts in high M stel i ar galaxies, and a lack of events 
with small A 2 blueshifts in low A/ 8te uar galaxies. This is in agree- 
ment with the Can results. Similarly to the Can H&K correlation, 
after correcting for the /i 2 -stretch correction, the significance of the 
relationship between the wavelength of A 2 and M ste uar drops to 0.5- 



4.7 Wavelengths of NUV A { and Si 11 4130 A features 

The position of the Aj feature (2920 A) was also measured. A 
weaker trend with stretch compared to A 2 was found, but in the 
same sense. No trends with B — V colour were found. The veloc- 
ity of the Sin 4130 A feature also correlates with stretch (higher 
stretch SNe have faster Sin velocities) at 2.3-<x. As discussed in 
Section |4~3l this feature is contaminated by contributions from S 11 
and Co n, and so may not be a clean tracer of the photospheric ve- 
locity. No trend between Sin 4130 A velocity and B - V colour is 
seen. We do not investigate the host galaxy nor Hubble residuals 
for these features since any trends are of lower significance, and 
instead focus our analysis on the Ca 11 H&K and A 2 features. 



4.8 Hubble residuals 

The Hubble residuals m'^" 1 (with B - V colour and stretch correc- 
tions) and m r g" > ' no " (with B- V colour, but no stretch correction) are 
shown in Fig.Q~3]as a function of the wavelength of the Can H&K 
feature for the HST sample. No trends between m" sid or m re B sii - n ° " 
and Can H&K velocity are identified in the data. This suggests 
that, at least for this small sample of 11 SNe, there is no addi- 
tional correlation between Ca 11 H&K velocity and Hubble residual. 
Fig.[l4]shows the Hubble residuals as a function of the A 2 position. 
Again, for m'™ d , no correlations are seen. However, a correlation 



with m^" d ' ,w 1 as a function of A 2 is apparent in the right panel of 
Fig-Elbut with low significance (1.7-cr). 



5 DISCUSSION 

In this paper, we have studied a large, low-z sample of SN la NUV 
spectra obtained with HST+STIS. This sample contains 32 high 
S/N NUV spectra around maximum light, greatly improving upon 
the size of previous samples. Photometric data were obtained for 
our sample using the LT and LCOGT robotic telescopes, so that 
light curve parameters could be calculated using the light-curve 
fitting routine, SiFTO. This allowed colour-correction of the HST 
spectra to account for the effects of dust extinction and intrinsic 
colour variations, a corr ection that could no t be applied in the ear- 
lier study of our dataset ( Cooke et al .120111) . 



5.1 Physical significance of an evolution with redshift and its 
dispersion 



In Section 13.21 we compared our new low-z mean spectrum 
to one constructed from the intermediate-z sample of E08, and 
found a modest, but statistically significant (3.1-cr effect), evolu- 
tion with redshift at NUV wavelengths (<3300 A), with the low-z 
mean spectrum having a depressed NUV flux compared with the 
intermediate-z mean. This was confirmed using measurements of 
the UV - b colour of the spectra in Section |4l4l where we showed 
that the mean UV — b colour of the low-z sample is 0. 1 8 mag redder 
than that of the intermed iate -z sample. Previous studies including 
E08, ICookeetalJ ( l201lh and lFolev et al] J2012al) also found a de- 
pressed NUV flux at low-z compared to intermediate-z. However, 
this is the first time that an unbiased, large low-z SN la sample has 
been used. We also find an increased spectral dispersion at NUV 
compared to optical wavelengths, confirming the results of previ- 
ous studies such as E08. 

Previous SN la modelling efforts have studied the ef- 
fect of metallicity and composition on SN la properties 
(e.g. Hoflich, Wheel er and Thiel emanrJ 1 19981; iLentz et alj 2000; 
ISauer et al] 2008). However, these studies have been unable to 
explain the size of the observed increase in dispersion at NUV 
compared to optical wavelengths in SNe la; they predicted much 
smaller variations than those observed in our NUV spectra. Some 
modelling efforts have focussed on asymmetries in the SN ejecta 
dKromer & Sim 2009), where viewing angle effects were investi- 
gated as a potential cause of this dispersion. However, these mod- 
els predict an effect of 0.5 mag in the V band, which is much larger 

than the scatter seen in optical studies of SNe la. 

However, the recent spectral synthesis models of Walker et al. 

42012b . which vary the metallicity of the SNe, shed light on this 
observed evolution in NUV properties with redshift. Using a se- 
quence of models at a fixed luminosity (log(Lb i/Lo) = 9.6), where 
the metal content is varied by a factor of 0.05-5 times relative to 
that of SN 2005cf, we integrate their model spectra through our UV 
and b filters, and find that the observed change with redshift in the 
UV — b colour of 0. 1 8 mag corresponds to a decrease in metallicity 
of 0.4 dex with increasing redshift. 

Galax y M slellar is a crude proxy for both gas-phase and stellar 
metallicity dTremonti et alj|2004l ; iGallazzi e t al. 20051) , Our sam- 
ples are matched in M ste uar ( see Section 12. IK however the re- 
lation between M stc n al . and metallicity is expected and observed 
to evolve with redshift (e.g. ISavaglio et alj|2005l ; lLamareille et al.1 
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Figure 13. The B - V colour and stretch-corrected Hubble residuals and colour-only corrected Hubble residuals for the HST sample as a function of Can 
H&K velocity is shown in the left and right panels respectively. The SN measurements are also split based on their stretch into low-stretch (red, open circles) 
and high-stretch (blue, solid circles). The dashed lines are the best-fitting lines to the data. 



l2009l : ICresci et al.l2012l) . To determine the expected metallicity de- 
crease with increasing redshift for our SN la samples, we firstly 
estimated the mean age of formation of the SNe at the different 
mean redshifts of the samples (z = and 0.6) using the star forma- 
tion history from | L i| 112008) and the t delay-time distribution from 
iMaoz et al] bOllh . These redshifts correspond to mean formation 
redshifts of 0.3 and 0.8 for the low- and intermediate-z samples, 
respectively. At these redshifts for fixed M ste iiar. we estimated a 
decrease in metallicity of ~0.2 dex from [Lamareille et al.l d2009h . 
To first order, this is not too different to the metal licity evolu- 
tion in ferred based on the spectral modelling study of Walk eret al.l 
d2012[) . 

These new SN la models also give insight into the origin of 
th e UV - b scatter. W e again compare to the UV SN la models 
of I Walkeret "all j2012l) to determine if variations in metallicity can 
explain this observed UV - b dispersion. We find that metallicity 
variation can result in a UV - b dispersion of 0.24 mag (larger 
than the observed 0.21 mag). Therefore, the increased dispersion at 
NUV compared to optical wavelengths is consistent with metallic- 
ity variations within the SN la sample, which manifest themselves 
at shorter wavelengths. 



5.2 Velocity shifts and the connection to SN la progenitors 

In Section[4] we found that 'normal' SNe la with broader than av- 
erage light curves (higher stretch) have NUV features with system- 
atically higher velocities than those with narrow light curves (lower 
stretch). This effect is most obvious in the Can H&K and NUV A 2 
(3180 A, labelled in Fig. [4j features, although a smaller effect is 
also present in the Sin 4130 A and A\ (2920 A) features. In Fig. 
1101 we have plotted the best-fitting line to our low-z data to show 
the correlation between Ca n H&K velocity and stretch, which has 
a significance of 3.4-cr. However, it is unclear if instead of a corre- 



lation, we are seeing a bi-modality in the data, with SNe with low 
stretches (s < 1.01) having lower Can H&K velocities and those 
with higher stretches (s > 1.01) having higher Can H&K veloci- 
ties. We independently confirm t he presence of these trends usin g 
the spectral fitting code, synapps dTh omas, Nugent &Mezall20 111) , 
where the model Can and Sin velocities outputted a re also found 
to show this trend of increasing velocity with stretch. lFolevN2012l) 
did not identify this trend between Can H&K velocity and stretch, 
probably due to the restrictive stretch cut of .? < 1.05 employed, 
which removes 7 high-stretch SNe la from their sample. 

Similar trends of increasing velocity dwells et al.l 1 19941; 
Fisher etal.1 1 19951 ; iMazzalil 1 19981 ; iMazzali et alj 120071) of NUV 
and optical spectral features with increasing light curve width (in- 
creasing luminosity) have also been identified in previous sam- 
ples of maximum light, as well as nebular phase, spectra. In 
these studies, the trend was observed when a much broader sam- 
ple of SNe la were studied, including sub-luminou s 91bg-like 
SNe la dFilippenko et al.lll992al;lLeibundgut et al.ll 19931) and over- 
luminous 91T-like SNe la i Filippenko et alj 1992bl) . Instead our 
sample focusses solely on a cosmologically useful SN la sample, 
and it is within this narrower sample that a trend of increasing NUV 
velocities with increasing light curve width (stretch) is observed. 

A possible cause of this correlation between velocity and light 
curve width could be differences in the kinetic energies of the SNe: 
high luminosity (high stretch) SNe la have larger 56 Ni masses, 
which may lead to larger kinetic energies, and hence higher ve- 
locities. Additionally, higher velocities lead to stronger line blan- 
keting in the NUV (due to more overlapping of the spectral lines), 
which could cause a depression of the UV pseudo-emissions, as is 
seen for our high stretch sub-sample in Fig. [5] High- velocity com- 
ponents (from detached shells of material) have been identified in 
most SNe la with spectr a obtained at phases earlier than one wee k 
before maximum light dMazzali et ai1l2007l; iBlondin et ai]|2012h . 
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Figure 14. The B — V colour and stretch-corrected Hubble residuals and colour-only corrected Hubble residuals for the HST sample are shown as a function 
of A 2 in the left and right panels respectively. The SN measurements are also split based on their stretch into low-stretch (red, open circles) and high-stretch 
(blue, solid circles). The dashed line in the right panel is the best linear fit to the data with a significance of 1.7-cr from zero. 



which may influence the measurements of spectral velocities at 
these early times. Blending of high-velocity features with photo- 
spheric features could result in higher overall velocities being mea- 
sured. However, the level of persistence of these high- velocity fea- 
tures to maximum light epochs is unclear and in the case of the 
Can H&K feature, can be confused with the Sin feature which is 
present in the same wavelength region. Therefore, we have used the 
synapps spectral models to distinguish between these scenarios and 
show that Can H&K velocity does correlate with stretch for our 
sample. 

Features at NUV wavelengths, including the A 2 feature, have 
been found using spectral synthesis modelling to be dominated by 
reverse fluorescence processes with the dominant species that con- 
tribute being Men , Sin, Sn, along wit h some Fe-group elements 
dSauer et al) l2008t IWalker et al.l |2012|) . The \ x and A 2 pseudo- 
emission features can be considered simply as regions of lower 
opacity, and will have bluer positions if the velocities and kinetic 
energies are higher. 

The mean spectra split by host galaxy M ste uar and also the 
analysis of the individual features show a trend of increased Can 
H&K velocity/blue-shifted A 2 feature for the low M 6te iiar galax- 
ies, or equivalently we see a lack of low-velocity/less blue-shifted 
events in low Mgteiiar galaxies, and lack of high-velocit y/mo r e blue - 
shifted events in high M 8te iiar galaxies, as claimed bv lFolevI d2012t) . 
However, we also show that this relation is driven by our identified 
trend between velocity and stretch, and the well-known correlation 
that brig hter (higher stretch) SNe la favour late-type (low mass) 
galaxies dHamuv et alJll99fj| , l200rj : ISullivan et alj200fj) . 

Higher velocity spectral features have also been found to be 
more common in SNe la that display outflowing material in high- 
resolutio n spectra, suggestive of a single- degenerate progenitor 
channel (Sternberg e t alJl20lTI : [Foley et al.ll2012bl) . This suggests 
that we can distinguish different progenitor scenarios using spectral 



feature velocities. In this paper, we find that higher stretch (higher 
luminosity) SNe la have, on average, higher velocities, and these 
SNe are more common in low M ste iiai galaxies. Following this ar- 
gument to its logical conclusion, we suggest that SNe la in low 
Msteiiar galaxies (predominantly young stellar populations), should 
result from the single-degenerate progenitor channel. Conversely, 
the lack of evidence of outflowing material in SNe la with low 
velocities and that form more commonly in high M stc u ar galaxies 
(predominantly older stellar populations), suggest that these pro- 
genitor systems may be linked to the double-degenerate progenitor 
channel. 

5.3 Cosmological implications 

The low-z mean spectra can be split by the value of the Hubble 
residuals, calculated from the light curve fits to investigate the ef- 
fect of NUV spectral variations on the use of SNe la for cosmology. 
We find no noticeable difference in the mean flux when the sample 
is split into positive and negative Hubble residuals. We also observe 
no correlation between the stretch- and colour-corrected, as well as 
only colour-corrected, residuals and the velocities/wavelengths of 
the NUV spectral features. This suggests that the observed NUV 
spectral variations do not directly translate into a cosmological bias. 
However, our observed NUV redshift evolution could have a sec- 
ondary impact on NUV-opt ical colours, whic h enter some distance 
estimation techniques (e.g. lConlev et alj|2008l) . 



6 CONCLUSIONS 

We have performed a detailed study of 32 near-UV (NUV) spec- 
tra of low-z (0.001 < z < 0.08) SNe la using HST+STIS. These 
spectra are complemented by gn-band light curves for the entire 
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sample from the LT and LCOGT robotic telescopes. We compared 
our low-z sample to that of the intermediate sample of E08 to test 
for evolution with redshift, as well as investigating correlations be- 
tween spectral, light curve and host galaxy properties within the 
low-z sample. Our principal conclusions are: 

(i) We observe an evolution with redshift of the near-UV con- 
tinuum at the 3-cr level, with low-z SNe la having depressed near- 
ultraviolet flux compared with the intermediate-z sample of E08 
(Fig.|4j- This is consistent with measurements of the UV-b colour, 
where the low-z sample has a 0.18 (3.6-<t) mag redder colour than 
the intermediate-z sample. 

(ii) We identify an increased dispersion at near-UV compared to 
optical wavelengths in our low-z SN la sample, as was found by 
E08 for an intermediate-z sample (Fig. [4). 

(iii) Using the SN la spectral synthesis models of Walk er et al.l 
d2012h . we have shown that this redshift evolution can be explained 
by a decrease in metallicity with increasing redshift of ~0.4 dex. 
The increased dispersion seen at near-UV compared to optical 
wavelengths can be explained using the same models by metallicity 
variations with the SN la sample. The observed redshift evolution 
is in agreement with galaxy evolution studies show a ~0.2 dex de- 
crease in metallicity over this redshift range. 

(iv) We find that SNe la with broader light curves (higher 
stretches) have systematically higher expansion velocities, in par- 
ticular seen in the Can H&K feature (Fig.[7J Fig. 1 10b . as well as 
blue-shifted NUV spectral features such as A 2 (Fig.l 12t . No corre- 
lation between velocity and B - V colour is found. 

(v) We show that SNe la in low M stc u ar host galaxies hav e 
high NUV spectral feature velocities, as claimed by Folevl J2012h . 
However, we attribute it to our velocity versus light curve width 
(stretch) correlation, and the well-known result that higher lumi- 
nosity, broader light curve SNe la are preferentially found in low 
M ste iiar galaxies (Fig.QjJ. 

(vi) We construct mean spectra of our low-z sample, split by 
Hubble residual, and find little variation in the mean spectra, indi- 
cating that NUV variations do not result in a clear SN la cosmology 
bias (Fig. [6). 

Our results underline the critical importance of the UV region 
for SN la studies. Given the qualitative success of the SN la models 
in reproducing the trends in our data, a more quantitative analysis 
using detailed spectra can be performed, to obtain important param- 
eters such as the density structure and elemental abundances of the 
individual SN events. This will enable us to put tighter constraints 
on the progenitor systems of SNe la, and how their progenitor chan- 
nels may vary with host galaxy properties. 



7 ACKNOWLEDGEMENTS 

MS acknowledges support from the Royal Society. AG-Y and MS 
acknowledge support from the Weizmann-UK "making connec- 
tions" programme. A.G. further acknowledges support by the ISF, 
BSF, a Minerva grant, the ARCHES award, and the Lord Sieff of 
Brompton Fund. DAH and BD are supported by the Las Cumbres 
Observatory Global Telescope Network. MMK acknowledges gen- 
erous support from the Hubble Fellowship and Carnegie-Princeton 
Fellowship. EOO is incumbent of the Arye Dissentshik career de- 
velopment chair and is grateful to support by a grant from the Israeli 
Ministry of Science. 

Based on observations made with the NASA/ESA Hubble 



Space Telescope, obtained at the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS 5-26555. 
These observations are associated with programs 1 1721 and 12298. 
The Liverpool Telescope is operated on the island of La Palma by 
Liverpool John Moores University in the Spanish Observatorio del 
Roque de los Muchachos of the Instituto de Astrofisica de Canarias 
with financial support from the UK Science and Technology Facil- 
ities Council. Observations were obtained with the Samuel Oschin 
Telescope at the Palomar Observatory as part of the Palomar Tran- 
sient factory project, a scientific collaboration between the Cali- 
fornia Institute of Technology, Columbia Unversity, La Cumbres 
Observatory, the Lawrence Berkeley National Laboratory, the Na- 
tional Energy Research Scientific Computing Center, the Univer- 
sity of Oxford, and the Weizmann Institute of Science. The William 
Herschel Telescope is operated on the island of La Palma by the 
Isaac Newton Group in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofsica de Canarias. Some of the 
data were obtained with the W. M. Keck Observatory, which is op- 
erated as a scientific partnership among the California Institute of 
Technology, the University of California and the National Aeronau- 
tics and Space Administration. These observations were made pos- 
sible by the generous financial support of the W. M. Keck Founda- 
tion. This paper uses observations obtained with facilities of the Las 
Cumbres Observatory Global Telescope. The Byrne Observatory 
at Sedgwick (BOS) is operated by the Las Cumbres Observatory 
Global Telescope Network and is located at the Sedgwick Reserve, 
a part of the University of California Natural Reserve System. This 
research has made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Laboratory, Cal- 
ifornia Institute of Technology, under contract with the National 
Aeronautics and Space Administration. 

This publication has been made possible by the participation 
of more than 10,000 volunteers in the Galaxy Zoo: Supernovae 
project, http : // supernova . galaxyzoo . org/ authors. 



REFERENCES 

Bailey S. et al., 2009, A&A, 500, 17 
Balland C. et al., 2009, A&A, 507, 85 
Benetti S. et al., 2005, ApJ, 623, 1011 
Blondin S. et al., 2006, AJ, 131,1648 
Bloom, J. S. et al., 2012, ApJ, 744, 17 

Blondin S., Mandel, K. S., Kirshner, R. P., 2011, A&A, 526, 81 
Blondin S. et al., 2012, AJ, 143, 126 

Brandt T. D., Tojeiro R., Aubourg E., Heavens A., Jimenez R., 

Strauss M. A., 2010, AJ, 140, 804 
Bronder T. J. et al., 2008, A&A, 477, 717 
Brown P. J. et al., 2010, ApJ, 721, 1608 
Bufano et al., 2009, ApJ, 700, 1456 

Cardelli J. A., Clayton G. C, Mathis J. S., 1989, ApJ, 345, 245 
Cenko S. B. et al. 2010, CBET, 2549 
Chomiuk L. et al., 2012, ApJ, 750, 164 
Chotard N. et al., 2011, A&A, 529, L4 

Conley A., Carlberg R. G., Guy J., Howell D. A. , Jha S., Riess A. 

G, Sullivan M„ 2007, AJ, 664, L13 
Conley A. et al., 2008, ApJ, 681, 482 
Conley A. et al., 201 1, ApJS, 192, 1 
Cooke J. et al., 2011, ApJ, 727, L35 

Cresci G., Mannucci F, Sommariva V, Maiolino R., Marconi A., 
Brusa M., 2012, MNRAS, 421, 262 



© 0000 RAS, MNRAS 000, 000-000 



NUV observations of SNe la 21 



D' Andrea C. B. et al., 201 1, ApJ, 743, 172 
Dilday B. et al., 2012, Science, accepted 
D'Odorico S. et al., 2006, SPIE, 6269, 626933 
Ellis R. S. et al., 2008, ApJ, 674, 51 (E08) 
Faber S. M. et al., 2004, SPIE, 4841, 1657 

Filippenko A. V., Richmond M. W., Branch D. et al., 1992, AJ, 
104, 1543 

Filippenko, A. V., Richmond, M. W., Matheson, T., et al. 1992b, 

ApJ, 384, L 15 
Fisher A. K., 2000, PhD thesis, University of Oklahoma 
Fisher A., Branch D., Hoflich P., Khokhlov A., 1995, ApJ, 447, 

L73 

Foley R. J. et al., 2008, ApJ, 684, 68 
Foley R. J. et al., 2008, ApJ, 686, 1 17 

Foley R. J., Sanders N. E., Kirshner R. P., 201 1, ApJ, 742, 89 
Foley R. J, 2012, ApJ, 748, 127 
Foley R. J. et al., 2012, AJ, 143, 113 
Foley R. J. et al., arXiv:1203.29T6V l 

Gallazzi A., Chariot S., Brinchmann J., White S. D. M., Tremonti 

C. A., 2005, MNRAS, 362, 41 
Gal- Yam A. et al., 201 1, ApJ, 736, 159 
Gonzalez-Gaitan S. et al., 201 1, ApJ, 727, 107 
Guy J. et al., 2007, A&A, 466, 1 1 
Guy L. et al., 2010, A&A, 523, 7 

Hachinger S„ Mazzali P. A., Benetti, S„ 2006, MNRAS, 270, 399 
Hamuy M., Phillips M. M., Schommer R. A., Suntzeff N. B., 

Maza J., Aviles R., AJ, 1996, 1 12, 6 
Hamuy M., Trager S. C, Pinto P. A., Phillips M. M., Schommer 

R. A., Ivanov V., Suntzeff N. B., 2000, AJ, 120, 1479 
Hoflich P., Wheeler J. C, Thielemann F. K., 1998, ApJ, 495, 617 
Hook I. M., J0rgensen I., Allington-Smith J. R., Davies R. L., 

Metcalfe N., Murowinski R. G., Crampton D., 2004, PASP, 1 16, 

425 

Horesh A. et al., 2012, ApJ, 746, 21 
Howell D. A. et al., 2005, ApJ, 634, 1 190 
Jin Z. et al., 201 1, CBET, 2708 

Kasen D., Ropke F. K., Woosley S. E., 2009, Nature, 460, 869 
Kelly P. L., Hicken M., Burke D. L., Mandel K. S., Kirshner R. P., 

2010, ApJ, 715, 743 
Kessler R. et al., 2009, ApJ, 185, 32 
Kramer M., Sim S., 2009, MNRAS, 398, 1809 
Lamareille F. et al., 2009, A&A, 495, 53 
Lampeitl H. et al., 201 1, ApJ, 722, 566 
Law N. M. et al. 2009, PASP, 121, 1395 
Le Borgne D., Rocca- Volmerange B., 2002, A&A, 386, 446 
Le Fevre O. et al, 2005, A&A, 439, 845 
Leibundgut B. et al., 1993, AJ, 105, 301 

Lentz E. J., Baron E., Branch D., Hauschildt P. H, Nugent P. E., 

2000, ApJ, 530, 966 
Li L.-X., 2008, MNRAS, 388, 1487 
Li W. et al., 201 1, Nature, 480, 384 
Lucy L. B., 1999, A&A, 345, 211 
Maguire K. et al., 2011, MNRAS, 418, 747 
Mannucci F, Delia Valle M., Panagia N, Cappellaro E., Cresci 

G., Maiolino R., Petrosian A., Turatto M., A&A, 433, 807 
Maoz D., Mannucci F, Li W., Filippenko A. V., Delia Valle M., 

Panagia N, 2011, MNRAS, 412, 1508 
Margutti R. et al., 2012, ApJ, 751, 134 

Markwardt C. B., 2008, in Astronomical Data Analysis Software 
and Systems XVIII, Non-Linear Least Squares Fitting In IDL 
withMPFIT,p251 



Mazzali P. A., Cappellaro E., Danziger I. J., Turatto M., Benetti 

S., ApJ, 1998, 499, L49 
Mazzali P., 2000, A&A, 363, 705 

Mazzali P. A., Ropke F. K., Benetti S., Hillebrandt W., 2007, Sci- 
ence, 315, 825 
Milne P. A. et al, 2010, ApJ, 721, 1627 
Nakano S. et al., 2011, CBET, 2783 
Nayak I., 2010, CBET, 2560 

Nugent P., Phillips M., Baron E., Branch D., Hauschildt P., 1995, 

ApJ, 455, L147 
Nugent P. E. et al., 201 1, Nature, 480, 344 
Ofek E. O. et al., 2012, PASP, 124, 62 
Oke J. B., Gunn J. E., 1982, PASP, 94, 586 
Oke J. B. et al., 1995, PASP, 107, 375 
Perlmutter S. et al, 1999, ApJ, 517, 565 
Pignata G. et al., 2009, CBET, 2022 
Rau A. et al., 2009, PASP, 121, 1334 
Riess A. G. et al., 1998, AJ, 116, 1009 
Riess A. G. et al., 2007, ApJ, 659, 98 
Russell B. R., Immler S., 2012, ApJ, 748, 29 
Sauer D. N. et al., 2008, MNRAS, 391, 1605 
Savaglio S. et al., 2005, ApJ, 635, 260 
Schaefer B. E., Pagnotta A., 2012, Nature 481, 164 
Schlegel D., Finkbeiner D., Davis, M., ApJ, 500, 525 
Silverman J. M., Ganeshalingam M., Li W., Filippenko A. V, 

2012, arXiv:1202:2130v2 
Silverman J. M., Kong J. J., Filippenko A. V, 2012, 

arXiv:1202:2129v2 
Smith A. et al., 201 1, MNRAS, 412, 1309 
Steele I., et al., 2004, SPIE, 5489, 679 
Sternberg A. et al., 2011, Science, 333, 856 
Sullivan M. et al., 2006, ApJ, 648, 868 

Sullivan M., Ellis R. S., Howell D. A., Riess A., Nugent P. E., 

Gal- Yam A., 2009, ApJ, 693, L76 
Sullivan M. et al., 2010, MNRAS, 406, 782 
Sullivan M. et al., 201 1, ApJ, 737, 102 

Timmes F. X., Brown E. F, Truran J. W., 2003, ApJ, 590, L83 

Thomas R. C, Nugent P. E., Meza J. C, 201 1, PASP, 123, 237 

Tremonti C. A. et al., 2004, ApJ, 613, 898 

Tripp R., Branch D., 1999, ApJ, 525, 209 

Tripp R., 1998, A&A, 331, 815 

Van Dokkum P. G, 201 1, PASP, 113, 1420 

Vernet J. et al. 201 1, A&A, 536, 105 

Walker E. S et al., 201 1, MNRAS, 410, 1262 

Walker E. S. et al., 2012, MNRAS, in press 

Wang X. et al., 2009, ApJ, 699, L139 

Wang X. et al., 2012, ApJ, 749, 126 

Wells L. A. et al., 1994, AJ, 108, 2233 

Yaron O., Gal- Yam A., 2012, |arXiv:1204.1891V 2 

York D. G. et al., 2000, AJ, 120, 1579 



APPENDIX A: WAVELENGTH AND VELOCITY 
MEASUREMENTS 



© 0000 RAS, MNRAS 000, 000-000 



22 K. Maguire et al. 



Table Al. Velocity measurements for Can H&K and Sin 4130 A, along with wavelength measurements for the Ai and Aj features. Both the measured values 
and the phase-corrected d measurements are given. 



SN name 


Measured 


Phase-corrected 


Measured 


Phase-corrected 


Measured 


Phase -corrected 


Measured 


Phase-corrected 




Can (10 3 kms~') 


Can (10 3 kms 4 ) 


Sin (10 3 kms~') 


Sin (10 3 kms~') 


X\ (A) 


A\ (A) 


An (A) 


Ao (A) 


PTF09dlc 


16.99±0.17 


17.78±0.67 


9.29±0.05 


9.52±0.17 


- 


- 


- 


- 


PTF09dnl 


16.81+0.15 


17.17±0.33 


10.09±0.03 


10.20±0.08 


2974±3 


2976±5 


3155.5±3.4 


3147.3±5.4 


PTF09dnp 


- 


- 


10.98±0.05 


11.51+0.35 


2952±3 


2962±16 


3168.1+3.3 


3132.2+15.6 


PTF09fox 


14.40±0.04 


15.14±0.62 


10.14±0.09 


10.38±0.18 


2963±20 


2967±21 


3160.2±4.8 


3143.9±8.6 


PTF09foz 


14.28+0.11 


15.07±0.67 


8.55±0.06 


8.81+0.18 


2976±5 


2982±9 


3170.3±4.6 


3152.8±8.9 


PTFlObjs 


16.47±0.22 


17.02±0.50 


12.69±0.06 


12.87±0.13 


2930±5 


2933±7 


3125.2±4.7 


3113.1+7.1 


PTFlOfps 


- 


- 


7.40±0.05 


8.37±0.62 


2927±4 


2946±29 


3181.9+5.3 


3117.2+28.1 


PTFlOhdv 


17.03±0.18 


17.94±0.78 


10.13 ±0.04 


10.42±0.20 


2962±5 


2967±1() 


3136.9+10.5 


3117.1±13.2 


PTFlOhmv 


15.02+0.11 


15.72±0.59 


9.29 ±0.05 


9.51+ 0.15 


2970±4 


2974±8 


3109.8±4.6 


3094.8±8.3 


PTFlOicb 


12.61+0.03 


12.85±0.20 


9.11 ±0.02 


9.19±0.05 


2953±3 


2954±4 


3154.2±2.3 


3148.8+3.1 


PTFlOmwb 


13.09±0.03 


12.99±0.09 


8.21 ±0.04 


8.18±0.05 


2949±5 


2948±5 


3174.4±3.2 


3176.2±3.5 


PTFlOndc 


14.20±0.12 


15.84+1.37 


10.18+0.15 


10.71±0.38 


2930±25 


2940±30 


3171.2±4.9 


3135.3+18.3 


PTFlOnlg 


15.20±0.17 


16.93±1.45 


9.32 ±0.09 


9.89±0.37 


2966±5 


2976±18 


3165.1±4.3 


3127.9±16.6 


PTFlOqjl 


14.88±0.19 


16.55±1.40 


9.33 ±0.06 


9.88±0.35 


2931+15 


2941 ±22 


3136.7±3.1 


3100.6+16.1 


PTFlOqjq 


11.87+0.08 


12.86±0.83 


10.02±0.05 


10.34±0.21 


2945±5 


2951+11 


3174.1±2.2 


3152.4±9.6 


PTFlOtce 


16.02±0.05 


17.00±0.81 


10.94±0.08 


11.25±0.22 


2931 ±3 


2937±10 


3160.6±5.3 


3139.7±10.4 


PTFlOwnm 


12.90±0.07 


14.04±0.96 


10.06±0.1 


10.442±0.27 


2947±37 


2954±12 


3164.1+2.1 


3139.0±10.9 


r Ir lUwot 




1 c /in_i_i i(\ 
1 J.49±l. 39 


n q ^ _i_o r\£. 
9.3j±U.U6 


n on_i_o 
9.89±U.35 


2989±5 


onnn_i_ 1 n 
lyyy+i 1 


31 /U.z±j.4 


3133. /±io.4 


PTFlOxyt 


15.01±0.08 


15.90±0.74 


10.86+0.1 


11.157+0.25 


2933±65 


2938±11 


3155.1+6.7 


3135.6±10.9 


PTFlOygu 


19.09±0.08 


18.99±0.12 






2928±10 


2927±10 


3155.3±3.9 


3157.1 ±4.1 


PTFlOyux 


16.05±0.13 


18.04±1.66 


6.97±0.14 


7.61±0.44 


2923±3 


2935±20 






PTFlOzdk 


16.95±0.14 


16.22±0.63 


8.30±0.04 


8.06±0.16 


2913±6 


2909±9 


3115.4+ 5.4 


3131.0±9.0 


PTFllkly 


11.93±0.04 


12.72±0.66 


8.54±0.13 


8.79±0.21 


2954±1 


2959±8 


3186.9+ 1.2 


3169.8±7.4 


SN2009le 


16.08±0.12 


16.16+0.14 


11.83+0.1 


11.86+0.11 


2948±31 


2948±3 


3148.2± 3.3 


3146.3+3.1 


SN2010ju 


16.95±0.07 


18.48+1.28 


9.99±0.01 


10.48±0.32 


2942±5 


2951+16 


3119.1+9.4 


3085.6±16.9 


SN2011by 


12.64±0.06 


12.56±0.09 


9.09±0.01 


9.06±0.02 


2958±2 


2958±2 


3166.5+1.1 


3167.9+1.3 


SN2011ek 


12.80±0.10 


13.85+0.88 


9.33±0.01 


9.67±0.22 
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